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Abstract 

We investigate the asymptotic behavior of the polynomials p, q, r of degrees n in 
type I Hermite-Pade approximation to the exponential function, defined by p{z)e~^ + 
q{z) + r{z)e^ — 0(2'^"+^) as 2 ^ 0. These polynomials are characterized by a Riemann- 
Hilbert problem for a 3 x 3 matrix valued function. We use the Deift-Zhou steepest 
descent method for Riemann-Hilbert problems to obtain strong uniform asymptotics for 
the scaled polynomials p{3nz), q(3nz), and r(3n2) in every domain in the complex plane. 
An important role is played by a three-sheeted Riemann surface and certain measures and 
functions derived from it. Our work complements recent results of Herbert Stahl. 
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1 Introduction 

The study of Hermite-Padc approximation for the exponential function was initiated by 
Hermite [19], in connection with his proof of the transcendency of e. A special case of this 
is the quadratic Hermite-Pade approximation of type I, where for any given three integers 
ni,n2,n3 > 0, one asks for polynomials Pni,n2,n3, <lni,n2,n3, and rni,n2,n3 of degrees ni,n2,n3, 
respectively, such that 

Pnun2,n3{z) + Qn^n^Me' + rni,n„n3 (^)e'' = z ^ 0. (1.1) 

These polynomials exist and are unique up to a common multiplicative factor. We put 



so that 
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Putting the left-hand side of (1.1) equal to zero, we obtain the quadratic algebraic Hermite- 
Pade approximants 



Qni,n2,n3{z) =1= '\/Qni,n2,n3{z)'^ ^j^ni,n2,n3 (■2^)^ni,n2,n3 (■2^) 
2rni,n2,n3(2) 

to e^. 

The analytic theory of the Hermite-Padc approximation to the exponential function was 
investigated, among others, by Mahler [26], Aptekarev [1], Chudnovsky [11], Borwein [10], 
Driver [16], Driver and Temme [17], and Wielonsky [41, 42]. The paper [27] by Nuttall deserves 
a special mention since it contains many deep ideas and conjectures about the asymptotics of 
Hermite-Pade approximants to general functions. Several surveys on the subject also exist, 
see De Bruin [9], Baker and Lubinsky [5], and Aptekarev and Stahl [2]. 

Very recently, Herbert Stahl [34, 35, 36] gave detailed asymptotic results for the scaled 
diagonal polynomials 

Pn{z) =Pn,n,n{^nz), Qn{z) = qn,n,n{^nz), Rn{z) = rn,n,n{^nz), (1.2) 

and for the remainder term 

En{z) = Pn{z)e-'^' + Qn{z) + Rn{z)e^''' (1.3) 

in the approximation. From these asymptotics, the limit distributions of the zeros of P„, Rn 
and En follow. The zeros accumulate on specific arcs in the complex plane. In Figure 1, the 
zeros of Pn,n,m Qn,n,n, and rn,n,n have been plotted for the value n = 60. It shows their re- 
markable distribution, especially the zeros of qn,n;n distribute themselves in a very particular 
way. These results on the distribution of the zeros of Hermitc-Pade approximants to expo- 
nentials may be seen as a continuation of the study initiated by Szego in [38] concerning the 
distribution of the zeros of Taylor sections of the series for e^, subsequently generalized by 
Saff and Varga in [31, 32, 33] to the zeros of the Pade approximants to e^, see also [22, 40]. 

Stahl [34, 35, 36] obtained his results from a saddle point analysis for the explicit integral 
formulas for Q^, and which are 

/c_i [w{w^ - 1)]"+^ 
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= ^li^-^-1)]-- ^'-'^ 

= [M^^- 1)1^+1 • ^'-'^ 

Here Cj is a closed contour in the complex plane encircling j in the positive direction, which 
does not encircle the other points in {—1, 0, 1}. The number C in (1.4)-(1.6) is a normalization 
constant. Stahl used C = n!2""'"^(3n)~" in order to make P„ a monic polynomial. For our 
analysis below, we found it convenient to make Q„ monic, and we choose 

C=^-^ . (1.7) 

(3ri)" ^ ^ 
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Figure 1: Zeros of the diagonal quadratic Hermite-Pade polynomials P60,60,60 (the diamonds 
on the left) (?60,60,60 (the circles in the middle), and r6o,60,60 (the boxes on the right) 



In this paper we propose a different approach to the asymptotic analysis. Our approach 
is based on a Riemann-Hilbert formulation for the polynomials P„, Qn, and combined 
with a steepest descent analysis for Riemann-Hilbert problems. This technique originated 
with Deift and Zhou [15] and was applied to the asymptotics of orthogonal polynomials by 
Deift et al. [12, 13, 14]. See also [7, 8, 4, 20, 21, 22, 23, 24] for recent developments. The 
Riemann-Hilbert problem for multiple orthogonal polynomials was given by Van Assche et 
al. [39]. In the present situation it gives rise to a Riemann-Hilbert problem for a 3 x 3-matrix 
valued function. As far as we are aware, this is the first time that the steepest descent method 
for Riemann-Hilbert problems is applied to a 3 x 3 problem. 

The Riemann-Hilbert problem is to find a 3 x 3 matrix valued function y : C \ F ^ C'^^^ 
where F is a closed contour in the complex plane encircling the origin once in the positive 
direction, such that 

1. y is analytic in C \ F. 

2. Y satisfies the jump condition 

/I ^-3n-2g-3n^ q\ 

F+(z) = F_(z) 1 , zeV, (1.8) 

\0 2-3n-2e3n2 i J 

where Y+{z) and Y-{z) denote the limiting values of Y{z') as z' approaches z G F from 
the inside and outside of F, respectively. 
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3. For large z 









-2n-2 



Z — ^ OO. 



(1.9) 








We will show in Section 5.1 that the Riemann-Hilbert problem has a unique solution and that 
i2i(^) = Pn{z), Y22{z) = z-3"-2Q„(z) for z outside T, Y22{z) = z'^''-'^ En{z) for z inside T, 
and 123(2;) =Rn{z). 

The steepest descent analysis consists of a number of transformations. To make the 
transformations work, we make heavy use of the works of Stahl [34, 35, 36]. A crucial role is 
played by the Riemann surface defined by 



which is considered as a three sheeted surface with cuts along two arcs Tp and Tp. The 
jumps of the different inverse mappings of (1.10) across the arcs determine probability mea- 
sures /Lip and supported on Tp and Tp, respectively. These measures turn out to be 
limiting distributions of the normalized zero counting measures of the polynomials Pn and 
Rn, respectively. 

We choose the closed contour T in the Riemann-Hilbert problem for Y so that it contains 
the cuts Tp and Tp. The measures hp and 11 p and their y-transforms 



are used to make the first transformation of the Riemann-Hilbert problem, which has the 
effect of normalizing the problem at infinity. Then we follow the general scheme, as presented 
in [14] or [12], for the asymptotic analysis of Riemann-Hilbert problems. It leads to a final 
Riemann-Hilbert problem whose solution has an explicit asymptotic behavior for n — > 00, see 
Theorem 6.4 in Section 6.3. Tracing our steps back to the original Riemann-Hilbert problem, 
we obtain strong asymptotic formulas for the Hermite-Pade approximants in every region of 
the complex plane, including the sets where the zeros are and their endpoints. 

In Section 2 we state the asymptotic results for the polynomials P„, Qn-, and and for 
the remainder as well as for the quadratic Hermite-Pade approximants. The asymptotic 
results make use of the functions obtained from the Riemann surface. The Riemann surface 
and the measures and functions derived from it are also described in Section 2. In Section 3 
we prove the statements about the Riemann surface and other geometrical objects involved in 
the problem. To prepare for the transformations of the Riemann-Hilbert problem we need a 
large number of relations between the various functions involved, such as the inverse mappings 
of (1.10) and the functions gp and gp. These properties are established in terms in Section 
4. Sections 5 and 6 contain the transformations of the Riemann-Hilbert problem and all 
asserted asymptotic results are finally proven in Section 7. 
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2 Statement of results 

2.1 The Riemann surface 

In order to state our results we first introduce an appropriate Riemann surface. The Riemann 
surface is chosen so that it parameterizes the critical points of the function 



w I— > 3zw — log [w{w'^ — 1)] . 
Note that the integrals in formulas (1.4)-(1.6) have the form 



i 



C, w{w^ - 1) 



„n ( 3zio -log [«; - 1 )] ) ^, 



W 



(2.1) 



(2.2) 



and that by the classical saddle point analysis for the asymptotic evaluation of integrals, the 
main contribution to the integral (2.2) comes from a critical point of (2.1). So we define TZ 
as the Riemann surface for the function 



z = z{w) 



1/1 1 1 

+ 7 + 



S \w w — 1 w + lj w{w^ — 1) 
Note that we obtain (2.3) if we set the derivative of (2.1) equal to zero and solve for z. 



np 



(2.3) 




Figure 2: The Riemann surface TZ 

The rational function (2.3) has three inverse mappings. These are the three solutions of 
the cubic equation 

1 



zw 



w 



ZW + - = 0. 



(2.4) 



The Riemann surface TZ consists of three sheets TZp, TZq, and TZr as shown in Figure 2, see 
Proposition 2.1. The bijective mapping i/; : TZ C is the inverse of (2.3). Wc denote its 
restriction to the three sheets by ipp, i/jq, and ipp, respectively. So tpp{z), iPq{z), and ipniz) 
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are the three solutions of (2.4). Typically wc will identify the three sheets with copies of the 
complex plane, and so ifjp, i/jq, and V'_R are defined on C with appropriate cuts. The sheets 
TZp and TZq are glued together along a cut Tp connecting two branch points zi and Z2, and 
the sheets TZq and TZp are glued together along a cut Tp connecting the other two branch 
points zs and Z4. 

The Riemann surface has four branch points zi = z{wi), Z2 = z{w2), z^ = z{w^)^ Z4 = 
z{wi), which are related to the points wi, W2, w^, ^4 for which z'{w) = 0. Simple calculations 
give 

Wk = 3-V4a;g-2fe-i, A; = 1, 2, 3, 4, (2.5) 

where the primitive 8th root of unity. The corresponding values of Zk = z{wk) 

are 

zi = 3-1/44^, Z2 = 3-^/^21 zs = 3-'/'uZ Z4 = 3-1/^4, (2.6) 
where ij02a = e2'^*/24 jg ^j^g primitive 24th root of unity. The precise sheet structure of TZ is 
given in the following proposition. 

Proposition 2.1. There is an analytic curve Tp from z\ to Z2 lying in the left half -plane, 
and an analytic curve Tp from Z3 to Z4 lying in the right half-plane, such that the following 
hold. 

(a) Three inverse mappings ipp, tpQ, and ipp of (2.3) exist so that ipp is defined and analytic 
on C\Tp, i/jQ is defined and analytic on C \ (Fp U Tp), and ipp is defined and analytic 
on C\Tp. 

(b) At infinity, we have the values ■0p(oo) = —1, iIjq{oo) = 0, and iljp{oo) = 1. 

(c) For z ETp, we have 

^,j^{^Q- i;p)+{s)ds e R, (2.7) 

with integration along the +side ofTp. [This is the side ofTp that is on the left while 
going from z\ to Z2 along Tp.] 

(d) For z G Tp, we have 

with integration along the -\-side ofTp. 

Because of symmetry, Tp is the mirror image of Fp under reflection with respect to the 
imaginary axis. We take the sheet structure on the Riemann surface TZ as in Figure 2. The 
functions tpp, xpQ, tpp are defined on the P, Q, and R sheet of TZ, respectively. Together they 
constitute a conformal map from TZ onto the Riemann sphere. The images of the different 
sheets are shown in Figure 3. 

Other curves of interest to our problem arc defined by the property that either ^ f^^ (^Q~ 
^pp){s)ds or ^ Iz^i^Q ~ '4'R){s)ds is real. These curves are described by the following propo- 
sition, see Figure 4. 

Proposition 2.2. (a) There are four analytic curves where ^ /^^(V'Q ~ ipp){s)ds is real. 
One of them is Tp, a second one joins Z2 with z\ and intersects the positive real axis. 
We call this curve Tp. The other two curves are unbounded and lie in the left half-plane. 
They join Z\ and Z2 with infinity, and we call them Te,i and Te,2, respectively. 



il)p)+{s)ds G 



(2.8) 
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56 of the Riemann surface Tl 





(b) There are four analytic curves where fz-^i^^Q ~ '4^r){^) ds is real. One of them is Tji, 
a second one joins with Z3 and intersects the negative real axis. We call this curve 
r|j. The other two curves are unbounded and lie in the right half-plane. They join 
and Z4 with infinity and we call them Te,3 and rE,4, respectively. 

(c) The curves T*p and Tr do not intersect. 

(d) The curves r|j and Tp do not intersect. 

(e) The curves T*p and r|j intersect in the points ±iy* on the imaginary axis. The value 
of y* is approximately y* = 0.621391 • • • . 

By symmetry, Tp and r|j are mirror images with respect to the imaginary axis. Similarly, 
Fg 1 and Fe4 are mirror images of each other, as well as Te,2 and Te,3- All contours are 
oriented as shown in Figure 4. The orientation induces a +side and a —side on each contour, 
where the +side is on the left and the —side on the right while traversing the contour according 
to its orientation. Propositions 2.1 and 2.2 arc proved in Section 3. 

Now we can also define the contours Fg and Fg. We define 

4 

= U TE,k. (2.9) 
k=l 

where Fe,^, k = 1, 2, 3, 4 are as in Proposition 2.2, and 

Tq = [-iy*, iy*] U (F^ n{Rez< 0}) U (F|j n {Re ^ > 0}) , (2.10) 

where Fp, F|j and iy* are as in Proposition 2.2. We choose an orientation on Tq by orienting 
the interval [—iy*, iy*] from —iy* to iy* and by orienting Fp from Z2 to zi and F^ from Z4 to 
Z3, see Figure 5. 
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Figure 4: Curves for which either J^^{iI)q - ^p)(s) ds or /^gCV'Q - '^r){s) ds is real. 



The contours Fp, Fq, Tr, and Fg divide the complex plane into six domains. We denote 
the unbounded domains by Doo,p, Doo,r, -Doo.c/ and -Doo,L) as shown in Figure 5. The bounded 
domains arc denoted by Dp and Dp, where Dp is in the left half-plane, and Dp is in the 
right half-plane, see also Figure 5. We also put 



-Doo = {Doo,P U Doo,R U D^^u U Doo,l) U (F^ \ {zi,Z2, Zs, Z4}) 



(2.11) 



This is the unbounded domain bounded by Fp, Fq and Fp. Of further interest is the bounded 
domain Dp bounded by Fp and Fp, see Figure 6, and its mirror image Dp with respect to 

the imaginary axis. 

This concludes the description of the Riemann surface. 

2.2 The measures fj,p, fj,Q, jip, and jiE 

We now define a measure on each of the curves Fp, Fq, Fp, and F^. The complex line 
element ds is defined according to the orientation of these curves given in Figure 5. 



Definition 2.3. We define a measure /ip on Fp by 



diJ,p{s) = -—{ipQ - tpp)+{s)ds s G Fp, 



(2.12) 



a measure /xp on Fp by 



dmis) = -2:^i'^Q ~ ds se Tp, 



(2.13) 
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r 1 1 r 

0.2 0.4 0.6 0. 



Figure 5: Curves Fp, Tq, Tr, and Tg, and domains Dp, Dr, D^^p,Doo,r,Doo,u, and Doo,L- 
Note that Tq does not only consist of the interval [—iy*,iy*] but also of the small arcs from 
Z2 via —iy* to and from 24 via iy* to zi. The domain Deo is the union of F^; with the 
unbounded domains Doqp, Doo,r, Doqu, and Doo,l- 



and a measure /xq on Fq by 



2m 



(V'Q - i^p){s)ds 



\ 27ri(V'K-'0P)(s)t^S 

The measure /ig on F^; is defined as 

3 

diiE{s) 



- i'R)is) ds 



s G r^n{Rez < 0}, 
s e F^n {Re^ > 0}, 
s G [-iy*,iy*]. 

s G Te,i U Fe,2, 
S G Te,3 U Te,4- 



(2.14) 



(2.15) 



A priori these are complex measures. Our first result is that fj,p, fxg, and hr are in fact 
probability measures and fiE is a positive measure. 

Theorem 2.4. We have that fip is a probability measure on Tp, jiq is a probability measure 
on Tq and hr is a probability measure on Tr. The measure /ie is a positive measure on F^. 

Theorem 2.4 follows from Propositions 4.2 and 4.8. 

The relevance of these measures is shown by the following theorem. For every polynomial 
p of exact degree n, we denote by Vp the normalized zero counting measure. Thus 



n f-^ 



p{z)=0 
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where each zero is counted according to its multiphcity. We also define a zero counting 
measure for the remainder function En, namely 

S„(z)=0 

where the normalization by n now corresponds to the degree of approximation and the 3n + 2 
interpolatory zeros of En at the origin have been excluded. 

Theorem 2.5. We have 

yPr, I^P, VQn I^Q, ^Rr. ^ (2-16) 

where the convergence is in the sense of weak* convergence of measures, i.e., jjin —>■ if 
J fdjin J fdii for every hounded continuous function f . Fkirthermore, we have 

z^E„ I^E, (2.17) 

in the sense that 

J f{s)di^E„{s) = J f{s)diJ,Eis) 

for every continuous function f such that f{s) = 0{s^^) as s oo. 

The convergence of the zero counting measures is due to Stahl [35]. The theorem shows 
that the measures hp, jiq, jiR, and he agree with the measures that Stahl introduced in a 
different way. 

In contrast to the measures fp^, i^q„, and up^ which are probability measures, the mea- 
sures have infinite mass. They also have unbounded support. As a result, the proof of 
the limit (2.17) is more involved than that of (2.16). 

2.3 The ^-functions 

For the strong asymptotic results we need the log-transforms (or complex logarithmic poten- 
tials) of the measures /xp, /xg, and fXR- 

Definition 2.6. We introduce three functions 

gp{z) = [ logiz -s)dnp{s), zeC\Tp, (2.18) 

JTp 

ggiz) = [ log{z-s)dfiQ{s), zeC\TQ, (2.19) 



9r{z) = / \og{z - s)diJ,R{s), zeC\rR, 



(2.20) 



which are defined modulo 2'Ki. 



Thus gp, gq, and gp are multivalued functions, depending on the specific choice of the 
branches of the logarithmic functions. Our results will involve expressions like e"^^ and e"^^, 
and then the multivaluedness will play no role. 
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2.4 Functions (fp and 



Two other important functions are the functions (pp and ifp given by 

3 




(2.21) 



and 



^Ri^) = o / ~ V'r)(s) ds 




(2.22) 



The paths of integration in (2.21) and (2.22) are in C \ (Fp U U {0}). The functions i^p 
and (pR are multivalued but the real parts are well-defined. Prom Proposition 2.2 we know 
that Recfp = on the curves Pp, Pp, P^',!, and rE,2, and that RepR = on Pp, P|j, Fe,!, 
and Te,2- We collect the main properties oi (pp and pn in the following lemma. 

Lemma 2.7. (a) T/ie rea/ part of pp is zero exactly on Pp, Pp, Pp,i, and Te,2- 

The real part of (pp is negative in D^^p U D*p, and it is positive in the remaining part 
of the plane. 

(b) The real part of pp is zero exactly on Pp, Pp, Pp,3, and Pp,4. 

The real part of pp is negative in -Dqo.r U Z)p, and it is positive in the remaining part 
of the plane. 

(c) On the imaginary axis, we have Repp = Repp. We have Repp < Repp in the left 
half-plane, and Repp > Repp in the right half-plane. 

Lemma 2.7 is proved in Section 4.5. 
2.5 Strong asymptotics away from the zeros 

Now we can state the strong asymptotic results for the polynomials Pn, Qn: and the 
remainder function En- Recall that the three polynomials are given by (1.4), (1.5) and (1.6) 
with the constant C as in (1.7) so that is a monic polynomial. Then P^i and Rn both 
have leading coefficients (— ^)"^^, and we will state our asymptotic results for the monic 
polynomials {—2)'^'^^Pn{z) and {—2)'^^^Rn{z). Throughout the rest of the paper, we will use 
the function \/3tt;^ + 1 which branches at the four points w^. given in (2.5). We choose as cuts 
for this function the two curves V'p+(rp) and ^p+(Pp) (see Figure 3), and assume that it is 
positive for large positive w. So, in particular we have that V 2>w^ + 1 = — 1 for = 0. 

The following theorem gives the strong asymptotics of the polynomials P„, Q„, Rn and 
the remainder term En away from their zeros. These results are due to Stahl [35], but we will 
give independent proofs below. 

Theorem 2.8. With the functions defined above, we have 




(2.23) 



uniformly for z in compact subsets of C\Tp 




(2.24) 
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0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 



Figure 6: The shaded region is where Re(pp is negative. It consists of the two parts -Dqo.p 
and Dp, where Dp is bounded by Tp and T*p. 



uniformly for z in compact subsets of C \ Fr, and 



Quiz) = < 



/ / 1 ■ 

= 1 + O - ) ) forzeD, 

' 1 + C ( - 1 1 for z^ Dp 



3ijj,{z) + 1 

^ngQ{z) 



l+O 



(2.25) 



for z e Dr, 



uniformly for z in compact subsets of C\ Fq . Furthermore we have 



1 \ " f,n{gR{z)+3z) / /I 



n 



for z G Doo,R 



'3i;%iz) + 1 

' 1 + O - ) ) forze D^,p 



3i/jUz) + 1 



En{z) = < 

^3V^(z) + l 
uniformly for z in compact subsets o/ C \ Fg. 



n 



^3n„-n{gp{z)+gR{z)) / /i 

1 + Q - 



/or z G Dp U U DocC/ U L>oo,L 
U Fq \ {2:1,2:2,2:3,^4}, 



(2.26) 
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2.6 Asymptotics near the curves Fp, Tq, Vr and Ve 

We will also obtain uniform asymptotics near the curves Fp, Fq, Fp and Fg, as well as in 
neighborhoods of the branch points z^. We start with the polynomials Qn and the remainder 
function En- 
Theorem 2.9. Uniformly for z in compact subsets of the region 

{Tq \ {zi,Z2, Z3, Z4}) U Dp U Dp U Doo,U U Doo,L, 

see Figure 5, we have 

,-2nipp{z) f fY 



Quiz) 



+ 1 



1 \ \ p-2n^R{z) 

i + o{-)) + 



Wr{z) + 1 



1 + C( - 

n 



. (2.27) 



Uniformly for z in compact subsets of Doo (which includes the arcs Te but not the branch 
points zi, Z2, Z3, Z4), we have 



En{z) 



= _2;^n^-n{gp{z)+gii{z)) 



+ 



-2n(pp (z) 



+ 1 



1 + 



n 



1 \ \ p-'^n<pji{z) 



Wniz) + 1 



1 + 



n 



.(2.28) 



Note that the asymptotic formula (2.27) (rcsp. (2.28)) holds in particular on Fq (resp. 
F^;), away from the branch points zi, Z2, 23 and 24, that is, on the curve where the zeros of 
Qn (resp. En) accumulate. 

It may be checked that the three terms in (2.27) are analytic continuations of the different 
asymptotic formulas wc have in (2.25). For z G Doo, wc have gqiz) = —gp{z) — g^lz) + 3log z 
by Corollary 4.11, and then (2.27) reduces to the first formula in (2.25) since Ileifp{z) > 
and Re(pR{z) > in Doo,u and Doo,l- For z G Dp, we have gQ{z) = -gp{z) - gR{z) + 
Slogz — 2ipp{z) and we obtain the second formula in (2.25). Finally, for z G Dp, we have 
9q{z) = —gp{z) — gR{z) + 31ogz — 2ipR{z), and wc get the third formula in (2.25). 

On Fq two of the terms in (2.27) have comparable absolute values. This causes the zeros 
of Qn to be close to Tq. At the points ±iy*, all three terms have comparable absolute values. 
Similarly, on F^ two of the terms in (2.28) have comparable absolute values, which causes 
the zeros of En in Dqo to be close to Fp;. 

For the polynomials P,,, and i?,„ wc have an asymptotic formula on Fp and Fp, respectively, 
away from the branch points, which now consists of two terms. 

Theorem 2.10. Uniformly for z in compact subsets of the region 



{rp\{zi,Z2})UD^,pUD*p 



14 



(this is the shaded region in Figure 6), we have 



+ 1 



l + O 



n 



2g2n(/?p(2) 



l + O 



n 



(2.29) 



where the +sign holds in DpU Dr and the —sign holds in D^o- 
Uniformly for z in compact subsets of the domain 



{Tr \ {zs, z^}) U {z G C I Re ^r{z) < 0} 



we have 



{-2r+^Rn{z) = e 



3V^(z) + 1 



1 \ \ 2p2"^«(^) 



S^*q{z) + 1 



l + 0[ - 

n 



(2.30) 



where the +sign holds in Dr U Dp and the —sign holds in 



For z away from Tp so that Re ^p{z) < 0, the asymptotic formula (2.29) reduces to (2.23). 
On Tp we have Reipp{z) = 0, and then the two terms in (2.29) are of comparable magnitude. 
For z G Fp, (2.29) can be re-written as 



{-2r+'Pn{z) 



Similar remarks hold for (2.30). 



2gnsp_ (z) 
SiPj.Jz) + 1 



1 + 



3V'4 + (z) + l 



l + O 



n 



2.7 Asymptotics near the branch points 

Near the branch points, the asymptotic formulas involve the Airy function Ai, which is the 
solution of the differential equation 



y"{z) = zy{z), 



satisfying 



Ai(.)^^-/..-|.-(l.o (_!,)) 



as z — !■ oo with | arg^| < tt, sec, e.g., [28] for more details. We only deal with the asymptotic 
behavior of P„, Qn and En near zi. Similar results can be given for the behavior near the 
other branch points. Here we will take the branch of the function <pp{z) which is at z = zi. 
So it behaves like 

^p{z) = c{z-z,f/^ + 0(^{z-z,f/^) 
as z ^ zi, with c 7^ 0. Then we define the function 



/i(^) 



2/3 



(2.31) 



15 



which is analytic for 2: in a neighborhood of zi. We take that 2/3rd power so that fi{z) is 
real and negative for z G Fp. Then 



fi{z) = ci{z-zi) + 0{{z-zif) 
as z ^ zi with some constant ci 7^ 0. Explicit calculations show that 

ci = f[{zi) = 2i/-^35/i2e-^". 
Theorem 2.11. There is a S > 0, such that we have, uniformly for \z — zi\ < 6, 



(2.32) 
(2.33) 



n^Hi{z) Ai ((n + 1)2/^1 (^)) (1 + ^ 

+n-V6^2(-^) Ai' ((n + lf/^fi{z)) U + O 



(2.34) 



(_2)'^+lQ„(z) = ^e^z^(n+l)(gp(z)+^p(z)-3z) 

nV6/,,(^)e-2-/3 Ai (e-2-/3(n + 1)^/^1 (^)) (l + C 

+n-V6/i2(^)e2-/3 Ai' (e-2-/3(„ + 1)'/Vi(^)) + ^ 

and 

i^/6/ii(z)e2-/3 Ai (^e2-/3(„ + 1)2/3/, (^)) (^1 + (^1^ j 
+n-V6/i2(^)e-2W3 Ai' (e2-/3(n + 1)^/^1 (^)) + ^ 



(2.35) 



(2.36) 



where hi and h2 are two analytic functions without zeros in\z — zi\ < 5, which have explicit 
expressions 

hi{z) = (iV2i(z) +iz-ie-3^iV22(z)) /i(z)i/4, 
with the branch of the fourth root in /i(z)^/^ taken with a cut along Tp, and 
h2{z) = {-N2i{z)+iz-\-^''N22{z)) fi{z)-^'\ 

Here 



2p-9pi.z) 

N21 iz) = , , N22 (z) 



+ 1 



V3V'^(z) 



+ 1 



(2.37) 
(2.38) 
(2.39) 



We use A^2i and A''22 to denote these functions, since in what follows they will appear as 
entries in a matrix N. Note that the function gp + (pp is analytic near z = zi. 

From the asymptotics near the branch points, one can deduce the behavior of the extreme 
zeros of Pn, Qn, Rn and near the branch points. We only state the result for the zeros of 
Pni Qn and En near zi. Recall that the Airy function Ai has only negative real zeros, which 
we denote by > — ii > —L2 > ■ ■ ■ > —l^, > ■ ■ ■ . 
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Corollary 2.12. Let z^.^, v = 1, . . . ,n, be the zeros of Pn, ordered by increasing distance to 
Z\ . Then for every u ^'H, we have 

= z, + 2-V33-5/i2e-i^i,^„-2/3 + o{^, (2.40) 

as n ^ oo. Let Zu,n, v = I, ■ ■ ■ ,n, be the zeros of Qn, ordered by increasing distance to zi. 
Then for every G N, 



zi - e'-^'l^^^.n-'^l^ + O 



zi + 2-V33-5/i2g-^-i,^„-2/3 + o(^^y (2.41) 



as n ^ oo. Let z^^^, v = 1, . . . ,n, be the zeros of En, ordered by increasing distance to z\. 
Then for every G N, 



z^ 

v,n 



f[{zi) \nj 

= Z, + 2-V33-5/12gi^i,^„-2/3 + O Q j , 



(2.42) 



as n ^ oo. 



2.8 Convergence of the Hermite— Pade approximants 

The polynomials P„, Qn, Rn satisfy 

Pn{z) + Qn(-2)e2"" + Rn{z)e^''' = O (^3^+2) , 

as z — >^ 0. Replacing the O term by 0, one might expect that the quadratic equation 

Pn{z) + Qn{z)Xn{z) + Rn{z)Xl{z) = (2.43) 

has a solution 

AKnyZ) 

which is close to the exponential e^^^ . This is indeed true for z in the domains Dp and Dr. 
Theorem 2.13. For z G Dp U Dr, there is a solution to (2.43) such that 

Xn{z) = e^^^(l + o(-\\. (2.45) 



The O-term is uniform on compact subsets of DpU Dr. 

Remark 2.14. The domain of convergence does not include the imaginary axis. On the 
imaginary axis one has roots of Qn — 'iPnRn which are branch points of accumulating 
along the interval [—iy*,iy*\. Hence, convergence does not take place in a neighborhood of 
even though the approximants come from a Hermite-Pade expansion around 0. 
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Remark 2.15. It turns out that dominates PnRn in the region Dp U Dp. If we agree to 
take the square root in (2.44) so that 

^/Qliz) - 4Pniz)Rn{z) ~ Qn{z), 

then in (2.44) we have to take the + sign in Dp and the — sign in Dp. So 



^ , , _ -Qn{z) + VQUz)-^Pn{z)Rn{z) 



Qn{z) + VQUz)-^Pn{z)Rn{z) Qn{z) 

and 



for z G Dp, 



^ -Qn{z)-y/QUz)-^Pn{z)Rn{z) ^ Qn{z) . , ^ „ 

Xn{z) = - p . s for z e Dr. 

Remark 2.16. It is possible to show that in D^c, there is no relative convergence to e^"^ of 
any of the two roots Xn{z). Absolute convergence holds in Dp and in a subdomain of Dp 
bounded by the curve Re {3z + 2(pp{z)) = 0. 

3 Geometry of the problem 

In this section we will prove Propositions 2.1 and 2.2 together. 
3.1 Proof of Propositions 2.1 and 2.2 

We recall that ipp, ipq, and ipp will be the three inverse mappings of 

z = Z{W) 



w(w^ — 1 



It is easy to see that z{w) is real if w is real, and that z{w) is purely imaginary if w G zR. In 
fact, for z G iR, there is one solution of the cubic equation (2.4) on the imaginary axis, one 
in the left half-plane, and one in the right half-plane. We denote these by i1^q{z), ijjp{z), and 
tpp{z), respectively. Then part (b) of Proposition 2.1 holds. We also note that ^q(O) = oo. 

These functions have analytic continuations, first to a vertical strip —e<Rez<e. around 
the imaginary axis. Then V'q has a pole at and from (2.3) it follows that 

tPQ{z) = ^ + 0{z) asz^O. (3.1) 



Near infinity, we have from analyzing (2.3) 



Mz) = + ^ + (3-2) 



1 „ / 1 



M^) = ^ + ^ 72 ' (3-3) 



Sz \z^ ' ^ 

3z~^ \ z 



Mz) = l + — + 0{^], (3.4) 
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as z —> oo. 



For < Rez < e, we then have that Jic^Q{z) > 0, Re^j:j(2;) > 0, while Re'0p(z) < 0. 
Since for z ^ iM, there can be no solution of (2.4) on the imaginary axis, it then follows by 
continuity that for every z in the right half-plane, there are two solutions of the cubic equation 
(2.4) with positive real part, and one solution with negative real part. For z with Rez > 0, 
we then define ipp{z) as the unique solution of (2.4) with negative real part. This defines ipp 
as an analytic function for Kez > —e. Similarly, for Rez < 0, we define i/jr{z) as the unique 
solution of (2.4) with positive real part, and we have V'R as an analytic function for Re2; < e. 

For the moment, we let Fp be an arbitrary simple curve that connects the branch points 
zi and Z2, and lies entirely in the left half-plane. We let Fp be the mirror image of Fp with 
respect to the imaginary axis and we assume that these cuts determine the sheet structure of 
the Riemann surface TZ. Then we can extend il)p, ipq, and V'i? to analytic functions on C\Fp, 
C \ (Fp U Fp), and C \ Fp, respectively. So we have part (a) of Proposition 2.1. 

Next we want to study curves so that ^ j^^{il)Q — il)p){s)ds is real. If z = z{t) is arclength 
parametrization of such a curve, then 



which upon differentiating leads to Re [z'{t){'4)q — 4>p){z{t))\ = 0. Since z'{t) / 0, and {%1)q — 
ipp){z) ^ (except at the branch points zi and Z2), we find that —z'{t)'^{'il^Q — ipp){z{t))'^ > 0. 
Thus the curve is what is known in geometric function theory as a trajectory of the quadratic 
differential 



see [29, 37]. Now, {tpQ —xpp)'^(z) is well-defined in the left half-plane, irrespective of the exact 
choice for Fp. It is analytic with simple zeros at the branch points zi and Z2, and a double 
pole at 0. 

Trajectories of the quadratic differential (3.5) which start from or end at zi or Z2 are 
called critical trajectories. From the local structure of trajectories of quadratic differentials, 
it is known that three trajectories emanate from any simple zero. So three critical trajectories 
emanate from zi, and we would like to label them by Fp, Fp, and Fp'^i. The labeling will 
depend on the global behavior of these trajectories. We consider the global behavior in the 
second quadrant, which we call G. 

It is known that any trajectory in G must begin and end either at zi, or at infinity, or on 
the boundary of G. It is also known that there can be no closed Jordan curve consisting of 
trajectories, since there are no poles of the quadratic differential in G. For these properties 
of trajectories of quadratic differentials, see [29, Chapter 8] and also [3]. So the three critical 
trajectories that start at zi, cannot end at zi, and they should end in G in one of the following 
ways 

(1) on the negative real axis, 

(2) on the positive imaginary axis, or 

(3) at infinity. 

If a critical trajectory meets the real axis, then by symmetry its mirror image with respect 
to the real axis will be its continuation to Z2- If there were two critical trajectories that come 



Re 




{4^Q-i:pf{z)dz\ 



(3.5) 
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to the negative real axis, then they both would continue to Z2- Then there would be a closed 
Jordan curve consisting of trajectories, which would not enclose a pole. This is impossible. 
Hence there can be at most one critical trajectory that meets the negative real axis. 

On the critical trajectories we have that ipp{z), see (2.21), is purely imaginary. Its deriva- 
tive is ip'p{z) = 2|?(V'Q - '^p){z), and 

Re<p'piz) = ^Imii^Q - i;p){z). (3.6) 

We claim that Iu^^|JQ{z) = ImV'p(-z) cannot happen for z on the positive imaginary axis. To 
prove this claim, let us assume that for z = iy with y > 0, we have Im^Q(z) = Iimpp(z). 
Note that ipQ maps the positive imaginary axis to the negative imaginary axis. So we have 
i/jQ{iy) = —iv with V > Q and 

v(v^ + 1) 

Recall that il^p{z) is always in the left half-plane. So we have tpp{w) = ~ ™ with u > 0, 
and by symmetry, we then have ■0ij(iy) = u—iv. Since ipp, ipg, and ipp are the three solutions 
of (2.4), we have 

ipp{z) + iPq{z) + tpR{z) = -, z e C. 

Then for z = iy, this leads to — 3w = i, so that y = However, this is incompatible with 
(3.7) ii V ^ 0. This contradiction proves the claim that ImipQ^z) = Imtjjp{z) docs not happen 
for z on the positive imaginary axis. As a result we see from (3.6) that the real part of (pp{z) 
is strictly monotone as z varies over the positive imaginary axis. Thus there can be at most 
one value y* > so that Reipp{iy*) = 0. Since Rec^p = on the critical trajectories, we see 
that there can be at most one critical trajectory that meets the positive imaginary axis. 

Now consider trajectories that end at infinity. At infinity, we have —{'ipQ — 'ipp)'^{(X)) = — 1. 
This means that the quadratic differential has a pole of order 4 at infinity, see [37], and all 
trajectories that extend to infinity arrive there with a vertical tangent. Assume that two 
critical trajectories extend to infinity in G. Then these trajectories are the boundary of a 
region G' in G. Any trajectory in G' begins and ends at infinity with a vertical tangent. 
However, according to [37, Theorem 7.4] a pole of order 4 has a neighborhood, so that any 
closed trajectory lying entirely in that neighborhood begins and ends at the pole, but from 
opposite directions. This is a contradiction, since there will be trajectories in G' that are 
arbitrarily close to infinity. This contradiction shows that there can be at most one critical 
trajectory that extends to infinity. 

So now we proved that each of the possible ways (1), (2), and (3) that a critical trajectory 
can end in G, happens at most once. Since there are three critical trajectories, it follows that 
each possibility happens exactly once. So there is a unique critical trajectory that meets the 
negative real axis. This trajectory continues to Z2, and we call it Fp. We then have that part 
(c) of Proposition 2.1 is satisfied. 

The critical trajectory that meets the positive imaginary axis is called Tp, and the critical 
trajectory that extends to infinity is called Te,i- By symmetry, we also find the trajectory 
rE,2 that emanates from Z2 and extends to infinity. 

By symmetry with respect to the imaginary axis, we also find the contours Tp, F^^s, rE,4, 
and in the right half-plane, so that for z on these contours we have that ^ Izsi''J^Q ~ 
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tpji){s)ds is real. This establishes part (d) of Proposition 2.1 and this proposition is now 
proved completely. 

Large parts of Proposition 2.2 have also been proved. All statements regarding Fp, Tji, 
^E,k, k = 1,2,3,4 are proved. We also know that Pp and P|j meet at the point iy* on the 
positive imaginary axis. What remains is to investigate how Pp and Pp continue beyond iy*. 
The following lemma will be useful. 

Lemma 3.1. The real part of (pp is strictly positive on Fr, Pe,3, and Te,4:- The real part of 

(fR is strictly positive on Pp, r_t;,i, and Te^2- 

We delay the proof of this lemma to Section 4.4, since its proof is easier once we have 
the functions gp and gu and know some of their properties. The reader may verify that the 
proofs that follow in Section 4 do not depend on the curves Pp and Pp. 

Assuming the lemma we can complete the proof of Proposition 2.2 as follows. Consider 
the quadratic differential —{ipQ — '4)p){z)'^dz^ in the domain G\ bounded by Pe,4, the curve 
Vr in the upper half-plane and the real interval (— oo,x*] where x* is the point where Pp 
intersects the real axis. The critical trajectory Pp has to end in G\ in one of the following 
ways 

(1) on the real axis, 

(2) on Pp or Ps,4, 

(3) at infinity. 

It cannot end on Pp or Pp 4 because of Lemma 3.1, and it cannot extend to infinity for 
the same reason that it could not extend to infinity in the left half-plane. So it has to end 
on the real axis, but it cannot end on the negative real axis, again for the same reason as 
before. So it has to end on the positive real axis, somewhere in the open interval (0, x*). 
Then by symmetry Pp continues in the lower half-plane to the branch point Z2 and Pp does 
not intersect with Pp. This proves part (c) of Proposition 2.2. 

By symmetry, also part (d) holds, and this completes the proof of Proposition 2.2. 

4 Measures and functions associated with the Riemann sur- 



In this section we study the measures /xp, //q, /ip, and /Ltg, and the functions (^p, (/?p, gp, gg, 
and gp. These measures and functions that are associated with the Riemann surface satisfy 
many relations that will be used in the transformations of the Riemann-Hilbert problem that 
follow in later sections. We also prove Theorem 2.4 and Lemmas 2.7 and 3.1 in this section. 

4.1 Properties of the measures //p, //p, and //g 
We start with a lemma. 
Lemma 4.1. We have 



face 




(4.1) 



and 




(4.2) 
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Proof. Let 7 be a closed contour on the sheet TZp going around Tp once in the positive 
direction. Then the residue theorem for the exterior of 7 gives 



because V'p is analytic outside 7 and we have (3.2). If we shrink 7 to Fp, then the integral 
becomes 

^.j^{{^p)^{s)-{^p)+{s))ds = l. 

Taking into account that (ipp)- = {'iI)q)+, we obtain (4.1). The reasoning is similar for the 
proof of (4.2), where we use a closed contour going around on TZr and the behavior (3.4) 
of ipR at infinity. □ 

Now we can prove the parts of Theorem 2.4 dealing with fip, and he- We give the 
proof of the part dealing with jjLQ only after we have a better understanding of the curve T*p 
and r|j, see Proposition 4.8. 

Proposition 4.2. /ip is a probability measure on Tp, ^xr is a probability measure on Tp and 

IJ,E is a positive measure on Te- 

Proof. The curve Fp is such that for z € Fp, the integral f^_^{ipQ — tpp)-{-{s) ds is real, see 
Proposition 2.1. For z = zi, it has the value 0, and for z = Z2 it has the value 1 by (4.1). Let 
z = z{t), t e [0,T], he arclength parametrization of Fp. The derivative of 



/ {iljQ-ijp)+{s)ds (4.3) 



2m 



is equal to ^{'4'Q{z{t)) ~ i^pizit))) z'{t) and this is different from for t G (0,T). Thus 
(4.3) is strictly increasing from for t = to 1 for t = T. This immediately implies that 
Hp defined by (2.12) is a probability measure on Fp. Similarly /xp defined in (2.13) is a 
probability measure on Fp. 

For He we observe that (2.15) defines a real measure on Fp, since, by Proposition 2.1, 
^Izi^'^Q ~ i^p)is)ds is real for z G Fp^i U Fp,2, and ^ J^.^iipQ - iJjR){s)ds is real for 
z £ Fp^3 U Fp;^4. Using an argument based on arclength parametrization, similar to the 
one above, we find that on each part Tej, j = 1,2,3,4, the measure fiE is either positive 
or negative. Since iPq{s) — ipp{s) = 1 + 0{l/s^) as s — 00 (see (3.2)-(3.3)), we have for 
z G Te,i, 

3 3 

{^Q -'>Pp)is)ds = —{z - zi) + 0{1) as z ^ 00, z eTE,i. 

Since Im{z — zi) — +00 as z — 00 along F^;^i, the integral is positive as 2; — 00 along 
Fp^i. As the measure ^up is of constant sign on Fp^i, we may thus deduce that it is positive 
everywhere on Fp^i. The reasoning for Fp^25 Fp_3, and Fp^4 is similar. 

This completes the proof of Proposition 4.2. □ 
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4.2 Properties of gp and gn 

The functions gp and gR were defined in (2.18) and (2.20). These are multi-valued functions, 
depending on the choice of the branch of the logarithm log(2; — s), which we assume depends 
on s € r in a continuous way. Since fj,p and fiR are probability measures, the gi-functions are 
defined modulo 27ri. 

Lemma 4.3. For the derivatives of the functions gp and gn we have 

g'p{z) = 3^p{z) + 3, zeC\rp, (4.4) 

and 

g'jiiz) = 32I;r{z) - 3, z G C \ r^. (4.5) 
Proof. The derivative of is easily obtained as 



3 f 1 

9'p{z) = 7r-- I (V'Q - '0p)+(s) ds. 

2m z — s 



If 7 is a closed contour going around Tp onTZp in the positive direction but with z outside 
7, then (since {'4>q)+ = {ipp)-) 



/ / X 3 / Vp(s) , 



The integral over 7 can be calculated with the residue theorem for the exterior of 7, for which 
there is a residue at z and at 00 given by (3.2). This proves (4.4). The proof of (4.5) is 
similar. □ 

Corollary 4.4. For z eC\{TpU Tr), we have 

25p(^) + 9'r{z) = ^ + 3 - S[Mz) - M^)], (4.6) 

and 

2g'Riz) + g'piz) = ^ - 3 - 3[^q(z) - Vi?(^)]. (4.7) 
Proof. We recall that ipp{z), iPq{z), and iPr{z) are roots of the cubic equation (2.4), so that 

i;p{z)+^Q{z) + i^R{z) = -^. (4.8) 

Hence by (4.4) and (4.5), 

'^9'p{z) + g'R{z) = 6il^p{z) + 6 + 3Mz)-3 

= 3[Mz) + + Mz)] - ^bpQiz) - Mz)] + 3 

= -^+3-3[iI;q{z)-Mz)]- 
A similar computation gives (4.7). □ 

A useful explicit expression of and gp in terms of the mapping functions tpp and ipp 
is given in the next lemma. 
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Lemma 4.5. The functions gp and qr have the following explicit expressions in terms of the 
mapping functions ipp and ipR, 



gp{z) = ^z{i,p{z) + 1) - \og[^p{z){^Pl{z) - 1)] - 1 + log - , /or z G C \ Tp, (4. 



and 



gniz) = 3z{iI;r{z) - 1) - logiM^Xi/jUz) - 1)] - 1 + log Q j , forzeC\ Tr. (4.10) 
Proof. We let 2; G C \ Fp and put w = '^p{z). Taking a derivative of (2.3), we find, since 



3 = - 



+ 



+ 



w"^ {w — 1)^ {w + 1)^ 



w'. 



Thus 



3^/jp{z) + 3 



-(w+l) 



+ 



+ 



u;-l)2 ('u; + l)2 



w 



1 1 

- + — + 

w 



+ 



+ 



w ~ 1 {w — 1)2 w+l 



w 



d_ 
dz 



logwiuP' - 1) -I h 



w w 



By (4.4), we then see that 



gp{z) = - log w{w'^ - 1) -I 1- 



w w 



+ C 



for some constant C. The constant can be determined from the behavior for z — >^ 00, since 
gp(z) = logz + C'(l/z), and w = il}p{z) = -1 + ^ + 0(1/^^). The result is that C = 2 + log ^ 
Using (2.3) we then find (4.9). 

The proof of (4.10) follows along similar lines. 



3- 

□ 



4.3 Jump properties of gp and qr 

We recall that the functions (fp and ipp were introduced in (2.21) and (2.22). The next 
lemma connects these functions with gp and gR. It will be frequently used in what follows. 
Throughout the rest of the paper we use ^ to denote the constant 



^ = log 2 — TTl. 

Lemma 4.6. For 2; G C \ (Fp U Tr), we have 

'^gp{z) + gpiz) = 3logz + 3z-2ipp{z)+e, 

'^■gRiz) + gp{z) = 3\ogz-3z-2^R{z)+L 
On the contours we have 

gp+{z) + gp-{z) + gpiz) = 3logz + 3z + i, zETr, 

9p{z) + gR+{z) + gR-{z) = 3\ogz-3z + t, z gTr, 



(4.11) 

(4.12) 
(4.13) 

(4.14) 
(4.15) 
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and 

gp+iz) - gp-{z) = -(pp+{z) + ifp-iz) = -2ipp+{z) = 2ipp_{z), z G Tp, (4.16) 
9r+{z) - 9r-{z) = -(PR+{z) + ipR-iz) = -2(pR+{z) = 2ipR_{z), z G r^j. (4.17) 
Proof. Integrating (4.6) from zi to z over some path in C \ (Fp U U {0}), we get 
29p(^) +9r{z) - 2gp{zi) - gR{zi) = Slog z - Slog zi + ?>{z - zi) - 2(pp{z), 
so that (4.12) holds with constant 

i = 2gp{zi) + gpizi) - Slogzi - Szi. 

Using the explicit expressions (4.9) and (4.10) for gp and gp we are able to show that i 
is equal to (4.11). To obtain (4.1S) we integrate (4.7) from Z3 to z. It turns out that the 
constant of integration is again given by (4.11). 
Next, we use Lemma 4.3 to find 

gp+{z) + gp-{z) + gR{z) = 3 (V'p+ + tpp- + ■iPr){s) ds + 3{z - zi), z G Fp. 

On Fp we have ipp-{s) = ^/^q+(s), so that 

gp+{z) + gp-{z) + gR{z) =3 {ipp + ipQ + V'ii)+(s) ds + 3z - 3zi, z e Fp. 

J 21 

Since V'p(s) + i^Q^s) + iPr{s) = ^, we obtain 

gp+{z) + gp-iz) +gR{z) = S log z - S log zi + 3z- 3zi = Slogz - 3z + £, z G Fp. 

This proves (4.14). A similar computation leads to (4.15). 

Finally, if we take (4.12) on the + and —sides of Fp and subtract, we get 

2gp+{z) - 2gp-{z) = -2^p+(z) + 2^pp-{z), z G Fp. 

This gives (4.16), since ippjf.{z) = —ipp-{z). Similarly we find (4.17). □ 

Since gp{z) = log 2; + 0{l/z) and gR{z) = log2; + 0{l/z) as z — 00, we get from (4.12)- 
(4. IS) that 

2(pp{z) = 3z + e + (-] asz^oo, (4.18) 



and 

2^r(z) = -3z + i + O (^^^ asz^oo. (4.19) 

We also see from (4.12) and (4.13) that (pp and ipp are multivalued functions, which are 
defined modulo ni, since gp and gp are defined modulo 27ri. 
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4.4 Proof of Lemma 3.1 



Now we are in a position to prove Lemma 3.1 which was needed in Section 3 to complete the 
discussion on the contours Fp and r|j. 

Proof. Let 7 be the contour which consists of Fji, Te,3 and Te,4- We define a function Gr 
on the right half-plane {Re z > 0} by 



This function is analytic across Tr because of the jump relation (4.17) of gR. The function 
u{z) = ReGR{z) is therefore harmonic in {Rez > 0} \ {Te,3 U Te,4)- On Te,3 and Te,^ we 
have ReipR = 0, and thus u is continuous in the full right half-plane. Furthermore KetpR < 
on the right of 7. Since gR is analytic across F^; 3 and F£;_4 it then easily follows that u is 
super harmonic in {Rez > 0}. Near infinity we have 



By symmetry Re gp{z) = Re gR{z) on the imaginary axis. Thus Re gp{z) — 3Re z is harmonic 
in the right half plane, Kegp{z) — 3Re z = u{z) on the imaginary axis, and u{z) — Kegp{z) + 
3Re2; is bounded from below as z — > 00. Then it follows from the minimum principle for 
superharmonic functions that u{z)—Iiegp{z) +3Re z > in the right half-plane. In particular 
this gives Re {gpi^z) — gp{z) + 3z) > to the left of 7 and Re {gR{z) — gp{z) + 2(pr{z)) > 
to the right of 7. Taking the difference between (4.13) and (4.12) gives gR{z) — gp{z) = 
—6z — 2ipr{z) + 2ipp{z). Then we see that Re (— 3z — 2ipR{z) + 2(pp{z)) > to the left of 7 
and Re {—3z + 2ipp{z)) > to the right of 7. Letting now 2^7 from either side, we get, 
since ReifR{z) = on 7, that 2Reifp{z) > SKez > on Tr, Fe^3, and Te,4- 

The proof that ReipR > on Fp, F^^i, and Fb_2 is similar, and also follows because of 
symmetry. □ 

4.5 Proof of Lemma 2.7 

Proof. By Proposition 2.1 and the definition of (pp, we have that Re(pp = on Fp, Fp, F^^i, 
and Fe,2- 

We know that Ret^p is a harmonic function in C\ (Fp U Fjj U {0}). Since iPq{z) ~ 1/z as 
z ^ 0; it easily follows from (2.21) that Iie(pp{z) — > —00 as z — 0. Then by the maximum 
principle for harmonic functions we get that Rc^^p < on Dp. 

As z ^ 00, we have (4.18). On the unbounded curves Te,! and Te,2 we have Re<y9p = 0. 
Prom (4.18) it then follows that these curves are asymptotic to the vertical line Re {Sz + i) = 
0, which is Rez = — ^log2. As z ^ 00 in the unbounded domain Doo,p we then have 
limsupRe(3z + ^) < 0, so that limsupRe99p(z) < by (4.18). Again it follows by the 
maximum principle for harmonic functions that Rec^p < on -Doo,P- 

For the remaining domain G = C \ {Dp U -Doo,p), we have that Reipp is harmonic on 
G\ Fp, with liminf (/9p(z) > as z — >■ 00 with z E G hy (4.18), and Re^^p > on Fp by 
Lemma 3.1. Thus again by the maximum principle, Re^pp > on G. This completes the 
proof of part (a) of Lemma 2.7. 




gpiz) to the left of 7, 

gR{z) + 2(pr{z) to the right of 7. 




Re logz -I- 0{l/z) as z — 00 to the left of 7, 

Re {-3z + log z + i) + 0{l/z) as z ^ 00 to the right of 7. 
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The proof of part (b) is similar and also follows because of symmetry. 

For part (c) we note that Re (fp = Re ipp on the imaginary axis because of symmetry. As 
z — >■ GO, we have by (4.18) and (4.19) that liminf Re ((^^^(z) — ipp{z)) > as z ^ oo in the 
left half-plane. Since Recpji and Re(pp are harmonic in {Re^ < 0} \ Fp and RecpR > Keipp 
on Fp, we find that Keipp < Keipp in the left half-plane. Similarly, we have Ret^p > Re(/?p 
in the right half-plane. 

This completes the proof of Lemma 2.7. □ 

4.6 The measure fj,Q 

Now that we have the full curves Fp and FJ^, we can prove the part of Theorem 2.4 that deals 
with HQ. We start with the analogue of Lemma 4.1. 



Lemma 4.7. We have 



and 



3 



3 
2^1 



[ {i,Q-i^p){s)ds = 2, (4.20) 
-J'^^{^Q-i^R){s)ds = 2. (4.21) 



Proof. If we deform the integral (4.20) over Fp to an integral over Fp, then we pick up a 
residue contribution at due to the pole of i/jq at the origin. So by (3.1) and (4.1) we have 



3 



T J^^ ii^Q - i'p)is) = ~ V'p)+(s) ds = 2. 



27rz ' JFp 

This gives (4.20) and in a similar way we obtain (4.21). □ 

Now we can prove that hq is a probability measure. 
Proposition 4.8. ij,q is a probability measure on Tq. 

Proof. We prove first that (2.14) defines a positive measure. On Fp we have that ^ J^^ {ipq — 
iljp){s)ds takes on real values, and has the value at z = 2:2 and the value 2 for z = zi by 
(4.20). Since its derivative does not vanish between Z2 and zi, wc find that ^{'ipQ — tpp){s)ds 
is a positive measure on Fp of total mass 2. Its restriction to Fp n {Re z < 0} is part of the 
measure //q, and so this part is a positive measure. Similarly, ^(ipR — il>Q){s)ds is a positive 
measure on F|j (due to (4.21)), and so also on Fp n {Rez > 0}. On the imaginary axis, we 
have, because of symmetry, that tpp{s) = —tpR{s), which means that tpR{s) — ippis) is real 
for s G iM. Then ^{il'R — tl'p){s)ds is a real measure on the imaginary axis. As s 00, 
we have '4>r{s) — ipp(s) = 2 + O(^). Since ipp — ipp never vanishes in C, this implies that 
^{i'R ~ il^p){s)ds is a positive measure on the imaginary axis, and so in particular on the 
interval [—iy*,iy*]. Hence all three parts of ij,q in (2.14) are positive, so that ij,q is a positive 
measure. 

Now we want to calculate the total mass of /ng . If we separate in J djiq the contributions 
from ipPjtpQ, and il^p, we get 
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Here stands for integration along the path from Zj to within Vq. We deform the 
contour from z\ to Z2 to the contour Tp. On Vp we have {ipp)+ = (V'q)-- Thus 



Similarly, we have 



Therefore, 



J ^^^^ '^Q{s)ds 

where 7 is the closed contour, consisting of the —sides of Fp and Tp and of the paths from 
Z2 to 2:3 and from Z4 to zi within Vq. Now we can use the residue theorem for the exterior 
domain of 7. The only contribution comes from infinity, and since iPq{s) = + ^{^) ^ 
s — > 00, see (3.3), we find Jp d/iQ = 1. Hence hq is a probability measure on Tq. □ 

4.7 The function qq 

For gq we have properties that are analogous to the properties of and gp obtained in 
Lemmas 4.3 and 4.5. 

Lemma 4.9. For the derivative of gq we have 



' 'Z) 



3ipq{z), z e Dc 

Sijp{z), z G Dp, (4.22) 

, 3'<Pr{z), z e Dp. 



Proof. We have 

9q{z) = I -^dfiq{s). 



z — s 



As we did in the proof of Proposition 4.8 to evaluate J dfiq{s), we can write the integral as 
an integral over the closed contour 7 that consists of the —sides of Fp and Fp and of the 
paths from Z2 to z^ and from Z4 to zi within Fq. We get for z G D^, 

Since ipq is analytic in D^o and vanishes at infinity, the only contribution comes from the 
pole at s = 2; and the result is that g'q{z) = 2>il)q{z) if 2 G L'oo- This gives (4.22) for z G D^. 
The expressions in the other two domains follow by analytic continuation. □ 



Lemma 4.10. gq has the explicit representation 

gq{z) = "izw - log[tt;('u;^ -!)]-! + log ^ '"^^^^ < 

Proof. This is analogous to the proof of Lemma 4.5. □ 



w = ^q{z) if z e Da 
w = ipp{z) ifzGDp, (4.23) 
w = 'iPr{z) if z e Dp. 
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We conclude this section with connections between gp, gg and gR that easily follow from 
the explicit expressions (4.9), (4.10), and (4.23). 

Corollary 4.11. (a) For z e D oo, we have 

9p{z) + gqiz) + gniz) = 31ogz. (4.24) 

(b) For z E Dp, we have 

9p{z) - gqiz) = 3z + 1 (4.25) 

(c) For z G Dp, we have 

gp{z) - gqiz) = -3z + i. (4.26) 
Proof. Adding the formulas (4.9), (4.10), and (4.23) with z G D^o, we get 

Sp{z) + gqiz) + gniz) = 

3z{wi +W2 + W3) - log [wi{wl - l)w2{wl - l)w3{w'i - 1)] - 3 + log (4-27) 

where wi = ipp(z), W2 = iPq(z), and W3 = iI)r{z). Since u^i, W2, and ^3 are the solutions of 
the cubic equation (2.4) we have 

z{w — w\){w — W2)(w — W3) = zw^ — VJ^ — zw + - 

o 

for every u> G C. Comparing coefficients of u)^, we get z{wi + W2 + w^) = 1, and taking 
It; = 0, 1/; = 1, and w = —1, we get ZW1W2W3 = — |, z{wi — l){w2 — l){ws — 1) = |, and 
z{wi + l){w2 + l)(ii'3 + 1) = |, respectively. Using this in (4.27), we arrive at (4.24). 

For parts (b) and (c), we start with the equation (4.24) in Doo- Making use of the jump 
relation (4.14) and (4.15) to go into Dp and Dp, we obtain (4.25) and (4.26). □ 

5 The Riemann— Hilbert problem and the first two transfor- 
mations 

Our asymptotic analysis is based on the Riemann-Hilbert problem for Y formulated in the 
introduction, see (1.8), (1.9). In this section we prove that the Riemann-Hilbert problem 
has a unique solution and that the solution is given in terms of the polynomials P„, Qn, and 
Rn, and the remainder En- We also do the first two transformation of the Riemann-Hilbert 
problem, which consist of a normalization of the problem at infinity, and a deformation of 
contours. 

5.1 The Riemann-Hilbert problem 

We show that the Riemann-Hilbert problem for Y has a solution in terms of the Hermite-Pade 
polynomials. 
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Theorem 5.1. Let Pn, Qn, Rn, o-nd En be as above. Then the solution of the Riemann- 
Hilbert problem for Y (see the introduction) is unique and is given by 



Y{z) 



Pn{z) 
\Pn,n~l,n+l{^nz) 



-3n-2 



9n+l,n-l,n(3n2;) rn+l,n-l,n{^'^z)\ 



-3n-2 



for z outside T, and 



Viz) 



/Pn+l,n^l,n{3nz) 

Pn{z) 
\Pn,n-l,n+\{'^nz) 



Z '^qn,n-l,n+l{^nz) 



z ^" 2e„+i,„_i,„(3nz) 



-3n-2 



Qn{z) Rn{z) 

rn,n-l,n+l{^nz) J 

rn+l,n-l,n(3nz)\ 

En{z) Rn{z) 



(5.1) 



,-3n-2 



(5.2) 



en,n-i,n+i(3n2;) rn,n-i,n+i{^nz) J 



for z inside T. In the first rows of (5.1) and (5.2) we use the Hermite-Pade polynomials of 
indices n+ l,n — l,n normalized so that Pn+i,n-i,n{'inz) is a monic polynomial, and in the 
third row we use the Hermite-Pade polynomials of indices n, n — 1, n + 1 normalized so that 
rn,n-i,n+i{^nz) is monic. 

Proof. The given Y is analytic inside and outside the contour T. This is clear from (5.1) and 
(5.2), except perhaps for the second column of (5.2) which has a possible singularity at the 



origin. However, the singularity is removable since e^^ 



0(^ni+n2+n3+2) as 2 ^ 0. 



The asymptotic condition (1.9) is satisfied because of the fact that p„+i,„_i,„(3nz), Qn{z), 
and rn,n-i,n+i{'inz) are monic polynomials. 

The jump condition can easily be checked. For the first and third entries in the first row 
it reads 

(Y-n)+ = (Yn)-, (^13)+ = (Yia)-, 
which is indeed so since Yn and Yis are both polynomials. We also have 



{yi2)+{z) 



z ^e„+i,„_i,„(3nz) 
^-3n-2 (p„^;^^„_;^^„(3nz)e" 
-3n— 2 —?>nz I 



+ gn+i,n-l,n(3nz) + rn+l,n-l,n{^nz)e^''^) 
3n—2„3nz 1 



z - ~e "--^(Fn)_(z) + (yi2)-(z) + z-"e^"^(yi3)-(z), 



and this is the jump condition (1.8) for the second entry in the first row. The second and 
third rows are handled in the same way. 

To prove uniqueness, we assume that Y is another solution of the Riemann-Hilbert prob- 
lem. First observe that det y is a scalar function which is analytic in C \ F. Because of (1.8), 
we have that (dety) + (z) = (dety)_(2;) for z G F, so that detl" has no jump, making detl" 
an entire function. For large z we have det 1^(2;) = 1 + 0(1^2;) by (1.9), hence by Liouville's 
theorem dety = 1 everywhere. We can therefore consider yy~^, which is analytic in C \ F. 
There is no jump on F since (yy~"^)+(z) = (yy^^)_(z) for every z € F, hence YY^^ is 
entire (i.e., each entry is an entire function). For large z we have YY^^{z^ = I + 0(1/ z), 
hence Liouville's theorem implies that YY^^{z) = I for every z, and hence Y{z) = Y(z). □ 

Remark 5.2. A Riemann-Hilbert characterization for general Hermite-Pade polynomials 
was given in 
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5.2 First transformation 

We will use the functions gp and qr, and the constant i = log 2 — m from Section 4.3 to 
transform the Riemann— Hilbert problem for y to a Riemann-Hilbert problem for U, given 

by 

/g-(n+l)ap(2) Q \ 

U{z) = L-^'-^Yiz) e^''+^')^<^P^^^+SR{z)] L"+\ (5.3) 

\ e-^n+l)gR{z) j 

where L is the constant diagonal matrix 

\ 

L = . (5.4) 

V e^/V 

For the contour F we take F = Fp U Fr U Tjj U F^, where Tjj is a contour connecting 
to zi and lying in Doo,Ui and F^ is a contour connecting Z2 to z^ and lying in -Doo,L- Then 
Rec^P > and Re^^R > on Tu U F^ by Lemma 2.7. 

We note that U is analytic on C \ F, since e^^^^-* and e^^^^-* are analytic and single- valued 
on C \ (Fp U Fp) and Fp U Fp c F. 

Since gp{z) = logz + C'(l/^) as z ^ oo, we have e("+i)ff^'(^) = z''+'^[l + 0{l/z)] as z ^ oo. 
Similarly we also have e("+^)^'«(^) = 2;"'+-'-[l + 0{l/z)] as 2; — oo. Hence 

C/(z) = / + oQj, z^oo. (5.5) 

So i7 is normalized at infinity. 

The jump relation for U needs to be worked out on the four pieces of the contour F = 
FpUFpUFc/UFi. 

For 2; G Fp we have 

(g- (n+l) [gp+ {z)-gp_ (2)] ^.gSz g(n+l) [-3 log z-Zz+gp+ {z)+gp_ {z)+gR (z)-E] q\ 
^{n+l)[gp+{z)-gp,{z)] q 

^^-^z^{n+l)[-i\ogz+iz+gp+{z)+2gR{z)-i] 1 

(5.6) 

Taking into account Lemma 4.6, we can simplify the jump (5.6) to 

(g2(n+l)</jp+(2) ^^Zz Q\ 

g2(n+l)¥'P-(2) I , Z G Fp. 

^g-32g-2(n+l)¥'fl(2) ij 

Similarly, the jump on Fp is 

U+{z)=U-{z){0 (.2{n+l)^R_{z) ), ^eFp. 

\0 ze-^^ 

On the parts Vu and F^, we have 

(1 2g32g(n+l)[-31og2-32+2ffp(2)+gfl(2)-^] qN, 
1 I . (5.7) 

2e"3^e("+^)["2'°s^+2^+^^(^)+2^«(^)"^l I J 
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If we use Lemma 4.6 then (5.7) can be re-written as 

1 01, z&Vu^Tl. 

^p-3z ^-2(n+l)tfR{z) 

Summarizing, we have the fohowing Riemann— Hilbert problem for U 

1. U is analytic on C \ F. 

2. U satisfies the following jump relations 

U+{z) = U-{z)\ g2(n+iVp,(2) zeTp, (5.8) 

^^-3z^-2(n+l)tpR{z) 
^e32g-2{n+l)<^p(2) 

U+{z) = U-{z) [ e2("+i)^«-W |, zeTR, (5.9) 

ze^'^^ 

_2g32g-2(n+l)¥'p(2) qN 

U+{z) = U-{z) \ 1 0|, zeTuUFL. (5.10) 

3. [/(z) = / + 0(i) asz^oo. 

In Figure 6 the shaded region is where Keipp is negative. For Reipu there is a similar 
picture, but reflected along the imaginary axis. The contours Fu and F^, are in the region 

where both Heipp and Re cpR are positive. The jump matrix in (5.10) for U on the contours Fjj 
and Fl is then the identity matrix / plus a matrix with entries that tend to zero exponentially 
fast as n — GO. 

The factor e-2("+i)¥'fl(2) in the last row of the jump matrix in (5.8) tends exponentially 
fast to because Re^?/? > on Fp, see Lemma 3.1. In a similar way the factor e~^("+^)'^^'^^) 
in the first row of the jump matrix in (5.9) on Fji tends to exponentially fast since Rc (/3p > 
on Fp (see Lemma 3.1 and Figure 6). Furthermore ipp+ = —(pp- is purely imaginary on Fp 
because of (2.7) and ipp^ = —ipp- is purely imaginary on Fp because of (2.8), so that the 
diagonal elements of the jump matrices on Fp and Fp are oscillatory. 

5.3 Deformation of contours 

The jump matrix in (5.8) can be written as a product of four matrices 

^g2(7l+l)v5p+(2) ^g32 qN^ 




ze^^ 



0> 

1 I . (5.11) 

2e~^^e~^("+^)'^«(^) 1, 
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Instead of jumping over Tp in one jump, we will make four smaller jumps, and rather than 
jumping over one contour, we jump over four contours, and each contour deals with one of 
the matrices in the product (5.11). We will open up a lens around Tp and introduce an 
extra contour for the last factor of the matrix. The lens consists of two contours Tp- (J rp+ 
connecting zi and Z2, such that Tp- is on the minus side of Fp and rp+ is on the plus side 
of Fp, but still inside the region where Heipji > and Re(/7p < 0. The contour rp++ stays 
away from zi and Z2- It starts at a point on the upper contour Tu and connects with a point 
on the lower contour Fl, so that it lies inside the region where Reipp > 0, (so it stays to the 
left of r^) , and does not intersect rp+ . The subarc of Tu between zi and the starting point 
of Tp++ we call Tjjp, and the subarc of from Z2 to the endpoint of Fp++ is called Tlp. 
These contours are drawn in Figure 7. 

1 



0.8- 

0.6 

0.4 

0.2 - 
0- 
-0.2 
-0.4 
-0.6 - 
-0.8- 

-1 




-1 -0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 
Figure 7: Deformation of contours around Fp. The dotted line is F|j 

In a similar way the jump (5.9) over Tr has a factorization 



1 



ze^^e 



2{n+l)<pp{z) 
g2(n+l)<pfl_ (2) 



ze 



3z 











1 







-l„3z 2{n+l)ipii+{z) 






-1 





-z ^e"^ 



_2g32g-2(n+l)ipp(z) 
1 





(5.12) 



We now open up contours Fp-, Fp+, and Fp++ around Fp, which are mirror images of the 
contours around Fp. The mirror images of Vyp and F^p are denoted by Vyr and F^r-, respec- 
tively. The remaining middle parts of Vy and F ^ are denoted by Fjjm and F , respectively. 
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-1 -0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 1 



Figure 8: Contours of the RHP for T 



All these contours are shown in Figure 8. All together there are 14 contours, and they 
determine 8 regions in the plane. The second transformation C/ i— > T will be defined in each 
of these regions separately. We define T as follows. 

We take 

T{z) = U{z), (5.13) 

for z in the unbounded region, and in the middle region bounded by Tp++, r^++, F^/m, and 
F^m. In the three regions near Fp we put 



Tiz) 



where 



U {z)Vp- (z), for z in the region bounded by Fp- and Fp, 

U{z)V-^+ {z)V-^ (z), for z in the region bounded by Fp and Fp+, 
^U{z)Vp^^{z), for z in the region bounded by Fp+, Fp++, Tjjp, and F^p, 

(5.14) 



1 0\ 

Vp-{z) = I ^-ie-32e2("+i)^^'(^) 1 , (5.15) 

1/ 

1 0\ 

Vp+{z) = I 2-ie-3^e2("+i)¥'pW 10, (5.16) 

1/ 

1 0\ 

Vp++{z) = I 1 . (5.17) 

,0 2e-^"^e-2("+^*'«(^) 1/ 
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In the regions near we put similarly 

T{z) 



U{z)Vji- (z), for z in the region bounded by F^- and Fr, 

U{z)V-^+ {z)V-^ {z), for z in the region bounded by Tji and F^+ , 
U{z)V~^^{z), for z in the region bounded by F^+, F^++, Tjjr, and F^r, 

(5.18) 



where 



'10 

Vr-{z) = I 1 2-^e3^e2("+i)^«W I , (5.19) 

^0 1 

'10 

Vr+{z) = I 1 z-'^e^z^2{n+l)^R{z) I ^ (5 20) 

^0 1 

'1 ze^^e~^("+^)'^-p('') 0^ 

Vr++{^) = I 1 I . (5.21) 

.0 F 



Then we have the following Riemann-Hilbert problem for T. 

1. r is analytic in each of the 8 regions, 

2. T has a jump on each of the 14 contours 

T+{z)=T_{z)Vs{z), zer. 



Si 



where s stand for any of the symbols P, P' , P+, P++, R, R- , R+, R++, U^, U"\ W, 
LP, L™, or U . The matrices Vp-, Vp+, Vp++, Vr-, Vr+, and Vr++ have already been 
defined above. The other jump matrices are 

/ ze^^ 0\ 
Vp{z) = -z-^e-^^ , (5.22) 
V 1/ 

/I \ 

Fij(z) = -z-^e^^ , (5.23) 
\0 / 

(1 2;e^^e~^("+^)'''^(^) 0\ 
1 , (5.24) 

1/ 

/I 0\ 

V[/r(z) = Fi.(z) = 1 , (5.25) 



Vu^{z) = Vl^{z) = \Q 1 I . (5.26) 

,0 2e-^^e-2("+^)'^«(^) 1> 



3. T{z) = I + 0{\) asz^oo. 
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Observe that all jumps, except for the jumps on Fp and Tr, tend to the identity matrix 
exponentially fast as n ^ oo. Hence we expect that the dominating contributions are the 
jumps Vp on Tp and Vr on Tji. 



6 Construction of parametrices and final transformation 
6.1 Pctrametrix for the exterior region 

We will now solve a Riemann-Hilbert problem for a matrix valued function A'^ on the contours 
Fp U Fp (see Figure 9) which, in view of what was said at the end of the previous section, is 
expected to describe the main contribution of the Riemann-Hilbert problem of T. 




1 1 1 r 

-0.6 -0.4 -0.2 



1 1 1 r 

0.2 0.4 0.6 0.8 



-1 -0.8 

Figure 9: Contours of the RHP for TV 

We look for : C \ (Fp U Fp) ^ C3x3 satisfying 

1. AT is analytic in C \ (Fp U Fp). 

2. N has jumps on Fp and Fp given by 



N+iz)=N^{z) 








-z 




zeFp, 



and 



'1 





N+{z) = N_{z) I 
,0 ze-^^ 








-z-^e^^ 



zeTp. 



3. N{z) = / + 0(i) as z^oo. 



(6.1) 



(6.2) 
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Proposition 6.1. A solution of the Riemann-Hilbert problem for N is given by 



/FiiM^)) F^iM^)) F,ii;Riz))\ 
N{z) = F2{^p{z)) F2{i^Q{z)) F2{iI^r{z)) I , 
\F^{^p{z)) F^{^q{z)) F-i{i^R{z)), 



(6.3) 



where 



Fi{w) = 
F2{w) = 
Fsiw) = 



—w{w — l)G{w) 
V3u;4 + 1 ' 
3(w^ - l)G{w) 

\/3wj4 + l~' 
w{w + l)G{w) 

V3w^ + 1~' 



(6.4) 
(6.5) 
(6.6) 



with \/3w^ + 1 defined and analytic in C \ {■ipp-^-{Tp) U iI^r-^-^Tr)), and such that it is positive 
for large positive w. The function G is defined by 



G{w) 



('m+l)(2'm-l) 
y;g w(w-l) 



for w G i/j{TZp) 



e w'^-i for w G ipClZq), 

(to-1)(2w+1) 



(6.7) 



I, we 



u(ui+l) 



for w G iljiUp). 



Proof. Let us consider the first row (A^ii, Ari2, -^13) of N. From (6.1) we get the following 
jumps on Tp 

r (A^ii)+(z) = -z-^e-^-{N^2)-{zl 

< (iVi2)+(^) = ze^'{Nri)-{z), zeVp, (6.8) 

[ (iVl3)+(^) = (iVl3)-(^), 

and from (6.2) the following jumps on Tr 



{Ni2)+{z) 



(A^ii)-(^), 

-3-(Ari3)_(z), 



ze 
—z 



Z^Tr. 



(6.9) 



-lg3z/ 



Clearly A'"!! is analytic on Tr and A^is is analytic on Fp. So we can (and do) see A^n as a 
function on the sheet TZp oi the Riemann surface TZ, N12 as a function on TZq and Ni^ as a 
function on T^-p. Then we transform the problem from TZ with the variable z, to the complex 
w-plane, via the mapping ijj : TZ ^ C The variables z and w are connected by (2.3). The 
images of the three sheets, the images of the branch points zi, Z2, Z3, Z4, and the images of 
the cuts Fp and Fp arc shown in Figure 3. 

Note that the images of Fp under the mappings V'P+ and tpp- (positive and negative 
boundary values of ■^p on Fp) give two arcs from wi to W2- They are oriented as shown in 
Figure 3. The orientation corresponds to the orientation of Fp. Together the arcs make up a 
simple closed loop around —1. 

Similar remarks hold for the two images ^^^+(1^^) and ^r-{Tr) of Fp. They make up a 
simple closed loop around 1. 
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Now we transplant the (as yet unknown) functions iVn, N12, and N13 from the Riemann 
surface to the w-plane, by defining Fi as fohows. 



Fiiw) 



Nu 











2 

w — 




w{w''' 













w G Tp{TZp), 



(6.10) 



Then Fi is analytic in C \ {ipp±(rp) U '0il±(rij)). The jumps that Fi should satisfy can be 
determined from (6.8)-(6.9) and are given by 



Fi+{w) = ze^^Fi_{w), 



Fi+{w) 
Fi+{w) 
Fi+{w) = ze-^^Fi_(u;) 



-z-^e-^^Fi_(w;), 
-z-ie^^Fi_(u;), 



w G V'p-(rp), 
w G V'P+(rp), 
u; G V;R+(rij), 
u; G iJjR-.{rR), 



(6.11) 



where z = z(i(;) = — r^-^- 

The asymptotic condition on A'' implies that iVii(z) — > 1, A^i2(^;) 
z — 00. For Fi, this means that 

Fi(-l) = l, Fi(0) = 0, i^i(l) = 0. 



0, iVi3(z) ^ as 



(6.12) 



We also want Fi{w) to have a finite limit as — ^ 00, since w = 00 corresponds to 2; = on 
the Q-sheet. 

We now seek Fi in the form 



Fiiw) = 



-w{w — l)G{'w) 



Then G should be analytic in C \ {ij;p±{rp) U iPr±{Tr)) with jumps 



(6.13) 



G+{w) = ze^^G-{w), 
G+{w) = z-hi-^'G-{w), 
G+{w) = z-^e^'G^{w), 
G+{w) =ze-^^G^{w), 

with z = z{w). The normalization for G is 

G(-l) = - 



w G V'p-(rp), 
w G V'p+(rp), 

w G tPr+{Tr), 
w G i/jR-iFp), 



(6.14) 



(6.15) 



It is straightforward to check that G given by (6.7) indeed satisfies (6.14) and (6.15). Then 
by (6.13) it follows that Fi has the correct jumps (6.11) and normalization (6.12). Then from 
(6.10) we recover A''!!, and Nis in terms of Fi by 

iVii(^) = Fi{^l;p{z)), Nu{z) = Fi{^q{z)), N,s{z) = Fi{^r{z)). 

Then the jumps (6.8) and (6.9) are satisfied, and in addition the normalization at infinity is 
correct. So we have found the first row of A''. 
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The proof for the second and third rows is similar. The only difference is that we have 
a different normalization at infinity, which leads to the construction of functions F2 and Fj, 
that satisfy the same jumps (6.11) as Fi, but are normalized by 



F2(-l) = 0, F2(0) = 1, ^2(1) = 0, 



and 



F3(-l) = 0, F3(0) = 0, ^3(1) = 1. 
Similar calculations then lead to the formulas (6.5) and (6.6) with the same function G. □ 

We remark that the entries of N have fourth root singularities at the branch points Zk 
(A; = 1, . . . , 4). More precisely, 

Lemma 6.2. The entries of N behave as follows near the branch points. As z ^ Zj with 
j = 1,2, we have 

(Nkiiz) =o{\z-z,\-y^), 



Nk2{z) =0{\z 

[iVfe3(z) =0(1) 

and as z ^ Zj with j = 3, 4, we have 



-1/4) ^ 



k = 1,2,3, 



(6.16) 



'Nk^{z) =0(1), 

Nu2{z) = {\z - Zj\-^/^) , 
NU^) =0{\z-z,\-y^), 



A; = 1,2,3. 



(6.17) 



Proof. Since w\ is a non-degenerate critical point of the mapping z = z{w), we have for the 
inverse w = ipp{z) as z — ^ zi = z{wi), 

i)p{z) =wi + c{z - zxfl'^ + 0{z- zi) (6.18) 

where c is a non-zero constant. Since wi is a simple root of 3w'^ + 1, it then follows that 



^J3iJj,{z) + l = C2iz - zxfl^ + O (|z - 



21, 



(6.19) 



with C2 / 0. Since the numerators of F\, F2, F3 as given by (6.4)-(6.6), do not vanish for 
w = wi, we find 

Nki{z) = FkiM^)) = 0{\z- z,\-^/^) , = 1, 2,3, 

as z ^ zi. In a similar way we find that Nk2{z) = O {^z — as z ^ zi. The entries in 

the third column of N arc analytic at zi, since tpfi is analytic at z\ and the functions F\, F2, 
F3 are analytic at iPr{zi). This proves (6.16) for j = 1. 

The behavior near the other branch points follows in a similar way. □ 

Remark 6.3. It will be useful to have another representation for the entries in the second 
row of N. They are 



N2i{z) 



2e~sp{z) 



+ 1 



7V22(Z) = - 



3V^4 (^) + 



iV23(^) 



+ 1 



(6.20) 
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It may be checked directly that these fTinctions have the right asymptotics as z ^ oo, and 
satisfy the correct jump relations on Tp and r^. They also satisfy the O-conditions of Lemma 
6.2. 

The Riemann-Hilbert problem for T is now very close to the Riemann-Hilbcrt problem 
for N because the jumps for T and A'' on the contours Tp and Tr are the same and the jumps 
for T on the other contours tend to the identity matrix as n — ^ oo, uniformly away from the 
branch points. So we expect that T behaves like N as n ^ oo away from the branch points. 
However, in order to justify this we need to analyze the Riemann-Hilbert problem in more 
detail near the branch points. 

6.2 Pcirametrices near the branch points 

We introduce four new contours Fi, T^, T4 which are small circles of radius S centered at 
the branch points, as indicated in Figure 10. We choose 6 small enough so that the circles 
do not intersect the curves Tp++ and r^++. The circles are oriented clockwise. Inside each 
of these contours we will solve the Riemann-Hilbert problem for T exactly. The analysis is 
similar for the four branch points, and we will only work out the analysis near zi in detail. 




-1 -0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 1 



Figure 10: Small circles Fi, F2, F3, F4 around the branch points 

Zooming in near the branch point zi gives a Riemann-Hilbert problem with five contours 
Fpj, Fp, Fp+, Ft/, Fi, as shown in Figure 11. 
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Figure 11: The contours near the branch point zi 
The jumps on these contours are 

/ 1 0^ 

Vp-{z) = Vp+{z) = L-ie-3^e2("+i)*'^W 1 

V 1; 

/ ze^^ 0> 
Vp = -z-^e-^^ 

V 1, 

(1 2:e^2e~^("+-^)'^-f 0^ 
1 

1, 

Wc look for a 3 X 3 matrix valued function M^^) defined within the disk Ai surrounded 

by Fi, such that 

1. M(i) is analytic in Ai \ (Fp U Fp- U Fp+ U Fj/p), 

2. M^^^ has the jumps 

M^^\z) = m'^}\z)Vs{z), zeVs, 
where s stands for any of the symbols P, P~ , P+, and U^. 

3. On Fi we have that M^^) matches N in the sense that 
uniformly for z G Fi. 
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The construction of M^^^ is basically a 2 x 2 problem, since the jump matrices arc non- 
trivial only in the 2x2 left upper block. The jumps have a standard form and a local 
parametrix can be built out of Airy functions. To be more precise, we will find M^^) in the 
form 

M«(z) =E«(z)^'W((n+l)2/Vl(2)) I ^l/2g3^/2g-(n+lVp(^) 

V 1, 

(6.21) 

Here we take the branch of fp which is zero at zi. So for z G Ai, we take 



3 

fp{z) = 2 / (V'Q - '<pp){s)ds 



with integration from zi to z in Ai \ Tp. Recall that ipp is defined modulo wi, and that so 
far the precise branch did not matter, since we only had expressions like e'^^'^'^^^'^^^^\ Since 

V'p(s) = wi + c{s - zxfl'^ + O (s - zi) , = wx- c(s - + C (s - 2l) , 

as s ^ zi, with the same non-zero constant ci, see also (6.18), we then have from this 
definition of </?p that 

^P{z) = {z- z{fl'^h{z), zeAi\ Fp, 

with h analytic and without zeros in Ai. The function (z — zi)^^"^ is defined with a branch 
cut along Fp. 

The function /i (z) is defined such that 

^Piz) = l[hiz)f\ 

It is a conformal map from Ai onto a convex neighborhood of 0. (We may have to shrink Ai, 
if necessary.) It maps Tp onto a part of the negative real axis. We still have some freedom in 
the choice of Fp- , F p+ and Tup . We take Tup so that /i maps Tup to a part of the positive 
real line. This means that Tup is an analytic continuation of Fp. Fp- and Fp+ are chosen so 
that they are mapped by fi onto rays in the complex s-plane. We denote the images of Fp, 
Fp-, Fp+, and Tup, by Sp, Ep-, Ep+ and S[/p. These contours are shown in Figure 12. On 
these contours we use the constant jump matrices 



Vp- = Vp+ 




on Tip- and Sp+, 



on Ep, and 



Vp 



Vup = 
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on Tijjp. This Riemann-Hilbert problem is well-known when we deal with 2x2 matrices (see, 
e.g., [12, §7.6], [13], [14, p. 1522]). Note however that the arrows are pointing in the opposite 
direction compared to Figures 7.18-7.20 in [12], so that we need to modify the solution a little. 
The solution ^'^-^^ uses the Airy function Ai{z), which is the unique solution of y"{z) = zy{z) 




Figure 12: The contours for ^''^^^ 
with asymptotics as z — oo given by 

which holds for | arg2;| < tt, where z^/^ and z^^^ are defined with principal branch (i.e., with 
a cut along the negative real axis). We will define ^^^^ by 

/Ai(s) -Ai(a;|s) 0^ 
*W(s) = Ai'(s) -a;|Ai'(a;|s) o|| e^'^/e o|, sel, 
V 1, 

/Ai(s) -Ai(w|s) 0\ /e-*^/6 0^ 
*W(s) = Ai'(s) -ulAi'iuls) 0| I e^'^/e o\Vp-, sell, 
\ 

(Ai(s) a)| Ai(w3s) 
Ai'(s) Ai'(a;3s) o|| e'""/^ o\v-^, selll, 


'Ai(s) Lv'^Ai^LUss) 
^'(^)(s) = I Ai'(s) Ai'(w3s) 0| I e*^/6 o|, s e IV, 


where (J3 = e^'^*/^ is a primitive third root of unity. 

With the above definitions of ^^^^ and /i it may then be shown that for any analytic 
prefactor E^^^ the matrix M^^^ defined by (6.21) satisfies the jump conditions on Tg, where s 
is any of the symbols P, P", P"*", and U^. The extra factor E^^^ has to be chosen in such a 
way that M^^^ satisfies the matching condition on Fi as well. 

On the part of Fi that lies in (/i)~^ (/) (the arc between Fp- and F[/p) the asymptotic 
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expansions for the Airy functions (6.22)-(6.23) give 
Ai((n + 1)2/3/1 (^)) ^ 



Ai((n + l)2/3a;|/i(z)) 

Ai'((n + 1)2/3/1 
Ai'((n + l)2/3a;i/i(z)) 



^(n + 1)-^/^ (/i(z))-^ e-("+i)-H^) (^1 + o (i)) , 

-^(n + l)-l/6 ihiz))'^ ^in/6^{n+l)<pp(z) 

2v7r V 



20F 
-1 



L(n + l)V6(/,(^))ie-("+i)^H^) (^l + o(i)) , 
^(n + l)V6 (_^^(^^)ig-m/6e(n+i)cpp(^) 1^1 + 1^1^^ 



Here the fourth root in (/i(2:))4 jg defined with a cut along Fp. This means that on this part 
of Fi we have 



mW(z) = 2^£;W(z) 



-in/6 



1)2/3/, (^))-i 0\ 

((n+ 1)2/3/1 (z))i 

1/ 







_g-i7r/6 _gi7r/3 q 
ly 









(6.24) 



In order to have the matching condition we therefore choose 



^^l/2g3^/2 ON 

z-V2e-3./2 
1, 

-in/6 _gW3 0\ "V((" + l)^^Vl(^))^ 







0' 



_g-i7r/6 _gi7r/3 g 

1> 




((n + 1)2/3/, (z))-i 

1, 

^l/2g3./2 0> 

z-V2e-3^/2 

1> 



((n+l)2/3/i(z))4 

((n + l)2/3/,(^))-i 

1/ 



(6.25) 



With this choice of E^-^^ it is then easy to see that the matching condition holds on the part 
of Fi in the region bounded by Tp- and Tjjp. A similar analysis for the other parts of Ti 
shows that the same E^^^ also works in the other regions. 

Prom (6.25) it is easy to see that E^^^ is analytic in Ai \ Fp. On Fp, both N and (/i)^ 

have a jump. N has the jump (6.1) and (/i)^ satisfies (/i)+ = —i{fi)t- Straightforward 

calculations then show that E^^ = E^^ on Fp, so that E^^^ is analytic across Fp. Prom (6.25) 
and the fact that the entries of N have at most fourth root singularities at zi, see (6.16), we 
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see that the entries of E^^^ have at most a square root singularity at zi. Since E^^^^ is analytic 
in Ai \ {^^i}, the singularity at zi is removable, and this proves that E^'^^ is analytic in the 
full Ai. This completes the construction of the parametrix M^^) in the neighborhood Ai of 

Zl. 

Prom the definitions (6.21) and (6.25) it is clear that 

(* * 0> 
* * 
1, 

where * denotes an unspecified entry. This means that the third column of M^^^ agrees with 
the third column of N. Wc will use this in what follows. 

In a similar way, we can construct parametrices near the other 

branch points Z2, Z3, Z4. The third column of M^^) agrees with the third column of N, and 
the first columns of M^^) and M^^^ agree with the first column of N. 

6.3 Third transformation 

We now introduce the final matrix 

r T{z) {N{z))-^ , z outside Ti, Fa, Tg, 

S{z) = I (6.26) 
[t{z) {M^\z)y^ , z inside Tj, j = 1,2,3,4. 

Inside each Tj the matrices T and M^^^ have the same jumps, hence S has no jumps inside 
Fj. Outside the Fj the matrices T and N have the same jump matrices on Fp and Fr. Hence 
S has no jump on Fp and Fji. This means that S solves a Riemann-Hilbert problem on the 
system of curves shown in Figure 13. 

S is analytic outside the above system of contours and it is normalized at infinity 

S{z)=I + o(-] , z^oo. (6.27) 



^z ^ 

Theorem 6.4. The matrix S{z) has the behavior 

S{z) = I + O (^^^ , (6.28) 

uniformly on C \ Sg, where S5 are the contours in Figure 13. 

Proof. The jumps on all of the contours are uniformly of the form / + 0(e^'^"') with some 
fixed c > 0, except for the jumps on the circles Fj where we have 

S+{z) = S-{z)M^^\z)N-^{z), z G Fj. 

Because of the matching condition we have 

M^^\z)N-^{z)=I + o{^ 

uniformly for z G Fj. Hence S{z) solves a Riemann-Hilbert problem, normalized at 00 
with jumps close to the identity matrix up to C(l/n), uniformly on the contours S5. We 
can then use arguments as those leading to Theorem 7.171 in [12] to obtain (6.28), see also 
[13, 14, 21]. □ 
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Figure 13: Contours of the RHP for S 



Using the asymptotic expansion for the Airy function and its derivative, it is possible to 
obtain a full asymptotic expansion for the jump matrix on the contours Tj in powers of 1/n. 
This in turn yields a full asymptotic expansion for S 



as in Theorem 7.10 of [14]. We will not use the full expansion, since (6.28) is enough for the 

asymptotic results as stated in Section 2. It will be clear however that the use of the full 
asymptotic expansion (6.29) will also lead to full asymptotic expansions for the polynomials 
Pn, Qn-, and Rn, and for the remainder En. 

7 Proofs of the asymptotic formulas 

We now know the asymptotic behavior (6.28) of S* as n ^ oo. We will trace back our steps 
to the original Riemann-Hilbert problem for Y to obtain asymptotics for the polynomials in 
the Hermite-Pade approximation problem for the exponential function. 

7.1 Proof of Theorem 2.8 

Proof. We start with the proof of the asymptotic formula (2.23) for P^. Let K be a compact 
subset of C \ Tp. We have the freedom to take the contours Fp-, rp+, rp++ near Fp, and 
the circles Fi and r2 around zi and Z2 in such a way that K is in the exterior of these curves. 
Let z & K. Then we follow the transformations Y ^ U ^ T ^ S. We see first from (5.3) 
that 




(6.29) 



Pn{z) = Y2i{z) = ?72i(z)e-("+i)^e("+i)^^(^). 
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Then from the definition of T in (5.13) and (5.18), we get that U2i{z) = T2i{z). We finally 
note that T{z) = S{z)N{z), if z is outside Ts and and T{z) = S{z)M(^\z) if z is inside 
Tj with j = 3,4. For T2i{z), we get in either case 

T2i{z) = S2i{z)Nn{z) + S22{z)N2l{z) + Si3{z)N3i{z) 

since the first columns of M^^^ and M^^^ agree with the first column of N. 
Since S = I + (i) , and since N21 does not vanish on K, we get 

P„(z) = iV2i(z)e-("+^)^e("+^)^^(^) (^1 + 

uniformly on K. Now we use the formula in (6.20) for N21 and we recall that i = log 2 — ni, 
to obtain (2.23). 

The proof of the asymptotic formula (2.24) for i?„ is similar. 

For Q n we proceed differently, because of the way Qn appears in the entries of Y . We take 
zo G C \ Fq, and show that there is a neighborhood A of zo such that (2.25) holds uniformly 
for 2: G A. First we assume that zq belongs to the outside region Doq. Then we can take the 
original contour F so that a neighborhood A of zq is in the outside region. Then for z e A, 
we have by (5.1) 

Q^{z) = z^''+^Y22{z) = U22{z)z^''+^e-^''+^^^3p{z)+9R{z))_ 

We can open the contours around Fp and Fji so that A is in the exterior region to these 
contours. Then we have U = T = SN={l + {^)) N, so that U22{z) = N22{z) (l + O (i)) 
uniformly for z G A. This leads to the first formula in (2.25), if we use (4.24) and the formula 
for N22 in (6.20). 

If Zo G Dp, then we can also open up the lenses around Fp and Tr so that a neighborhood 
A of Zq is not contained in these lenses. Then we have for z G A, by (5.2) 

Qn{z) = En{z) - P„(z)e-3- - R^{z)e^^' 

and we need to find out, what is the dominant contribution as n gets large. We already have 
asymptotic formulas for P„ and from which it follows that 

^ log |P„(z)e-=^"^| ^ Re {gp{z) - 3z - 1), (7.1) 

i log |i?n(-2)e=^"^| ^ Re {gR{z) + 3z - 1). (7.2) 
For En it is easy to obtain in a similar way 

-log|£;„(z)| ^Re(31ogz-5p(z)-5R(z)). (7.3) 
n 

Now it turns out that for z G Dp the term P„(z)e~^"^ dominates. Indeed, we have by (4.12) 
and (4.13) 

{gp{z) -3z-e)- igniz) + 3z-i) = {gp -gn- 6z) = 2{-^p{z) + ^r{z)) 
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and we know that the real part of —ipp + ^pR is positive in Dp. We also have by (4.12) 

{9p{z) -?>z- l)~ (3 log z - gp{z) - gR{z)) = 2gp(z) + gniz) - 31ogz - 3z - ^ = -2lpp{z) 

and the real part of —2ipp is positive in Dp. In view of (7.1)-(7.3) we then obtain the formula 
for Qn 

Qn{z) = -Pn{z)e-^''' {l + 0(- 



uniformly for z G A, and from what we already know about P„, 



(_l)n+l ngpiz) / /i 



n 



uniformly for z € A. This leads to the second line in formula (2.25) because gp{z) — 3z + 
log(— ^) = gqiz) for z G Dp, see (4.9) and (4.23). The proof for the case zq G Dp is similar. 

We finally have to consider the case that zq is on Fp or on Tp (but not one of the branch 
points). We consider only zq G Fp, since zq G Vp will follow in the same way. So let 
zo G Fp \ {zi,Z2}. After opening up lenses around Fp we have the contour Fp- to the left 
of zq, and the contour Fp+ to the right. We can take a neighborhood A of zq that is strictly 
contained in the domain bounded by F p- and Fp+ . Then for z G A H D^ , we have 

Q„(z) = z3"+2y22(^) = C/22(z)z3"+2e-("+i)(fl^'(^)+s«(^)) 

and U{z) = T{z)Vp}{z). Now we have 1/22(2) = ^22(2;), since the second column of Vp} 

is simply (O 1 O)'^. We open the circles Fi and F2 around zi and Z2 so that A is in the 
exterior. Then T = SN = {I + 0{^))N, and so 

Qniz) = ^3n+2^2^(^^^3n+2g-(n+l)(3p(.)+<;K(.)) (^J + Q (^^^^ 

uniformly for z G A fl D^. Using (4.24) and the formula (6.20) for Ar22, we then find (2.25). 
For 2; G A n Dp, we have 

Qniz) = Eniz) - Pniz)e-^''' - Rniz)e^^' 

= z^^+^Y22iz) - e-=^"^y2i^ - e^'^'Yniz) 

^3n+2^-{n+l){gp(z)+gR{z))jj^^(^^^ _ ^-3nz^{n+l)igp{z)-e)JJ^_^(^^^^ 
_ g3n2g(n+l) (Sfl {z)-e} U^^iz). 



= z 



Now we have that 
with 



Uiz)=Tiz)Vp+iz)Vp++iz) 



1 

Vp+iz)Vp++iz) = I 2-ie--^^e2("+i)^^'(^) 1 

^^-3z^-2{n+l)(pR{z) I 



48 



Then 



+ 



722 (^) 



+ 



^3n+2^-{n+l){gp{z)+gR{z))^^-3z^-2{n+l)v>Riz) _ ^3nz ^{n+l)igR{z)-i) 



Tisiz). 



The funny thing is that the expressions multiplying r22 (z) and Tis {z) are exactly zero. This 
follows from (4.12) and (4.13). Thus only the first term remains. This gives 

Q^{z) = -e9^W-^e'*^«(^)r2i(^) 

= -e»^(^)-^iV2i(z)e"^«(^) (^l + O (^i^^ 

= , ^ e"gg(^) (l + o(-)) 

uniformly for 2; G A PI Dp. In the last step we used (4.9) and the expression for N21. This 
proves the asymptotic formula ior z E An Dp. 

The asymptotic formula (2.26) for En is established in DpUDrUTq again by following the 
transformations Y U T 1— > 5. The asymptotics in the other domains may be obtained 
by observing that the main contribution in the sum 



,3nz 



is given by Pn{z)e Qn{z) or Rn{z)e^"-^ if z belongs to the domains Dao,p, D^^u U Doo,L 
or -Doo,R respectively. This completes the proof of Theorem 2.8. □ 



7.2 Proof of Theorem 2.5 

Proof. The limits for the counting measures vp^, uq^, and up^ follow from the strong conver- 
gence results established in Theorem 2.8. The proof using the unicity theorem for logarithmic 
potentials (see e.g. [30, Theorem II. 2.1]) is the same as the proof of [35, Theorem 2.1]. 

The proof for the limit of the measures is more difficult, since these measures have 
unbounded support and infinite mass. Our proof is different from the one given by Stahl [35]. 

Since En is an entire function of order 1 with a zero of multiplicity 3n + 2 at the origin, 
we have by the Hadamard factorization theorem 

En{z) = ^3n+2gA„2+B„ JJ / ^ _ ^ j gz/z,,„^ 

with Aji,Bn G C and where j = 1,2, .. . are the zeros of E^ different from the origin. 
The zeros are counted according to their multiplicities, and we assume that they are arranged 
in non-decreasing absolute value, so that < |2;i,n| < 1-2^2,71! ^ " " " ■ Taking logarithms, we get 

- log (z-3«-2e,(z)) = -{Anz + Bn)+ [ (log (1 - -) + -) di^E^s)- 



49 



Prom the strong asymptotic formula (2.26), we obtain from this that 



Hm 

n— +00 



^{AnZ + Bn)+ J (log (1 - ^) + J) dl^EAs) 



9e{z) (7.4) 



where 

gR{z) + 3z - 3 log z - ^ for z G D^^r, 



9e{z) = { 5rp(z) -3z-31ogz-£ forzGL>oo,P, (7.5) 



-9p{z) - QRiz) for z £ C\{Doo,rU Doo,p). 

Note that gE is analytic across Tp and Tr because of the relations (4.14)-(4.15), so that 
is analytic in C \ Te- The limit (7.4) holds uniformly for z in compact subsets of C \ Te- 
We first show that 

9e{z) = J (log (1 - J) + J) dMs) - 9P{0) - 9r{0) (7-6) 



for 2; G C \ r^;. To prove (7.6) we assume z G C \ {Dao,R U -Doo,p)- We take the derivative of 
the right-hand side of (7.6), which is 

The integral over F^;^! U rE,2 can be deformed to the integral over the —side of Tp, which 
leads to 



3 

271^ 



JrE,iurE,2\^- ^ 2TnJrp\z-s sj 

- <p ( — 1- - I ipp(s)ds 

ij^\z-s sJ 



_3_ 
2Tri 



where 7 is a closed contour going around Fp in the negative direction and with z and outside 

7, sec also the proof of Lemma 4.3. We compute the integral using the residue theorem for 
the exterior of 7, for which there is a contribution from the poles at z and 0. The result is 
that the first term on the right-hand side of (7.7) is equal to —3^pp{z) + 3V'p(0). In the same 
way, we find that the second term is equal to —3ijjR{z) + SiPr{0). Using (4.4), (4.5) and the 
fact that "0^(0) = y^1/3 = — 0p(O), we see that the derivative of the right-hand side of (7.6) 



is equal to —g'p{z) — Qr^z) for z G C \ [DooM U Doo,p)- Because of the definition of qe in (7.5) 
we then get (7.6) in case 2; G C \ {D^^r U D^^r). The formula (7.6) in the other domains 
follows by analytical continuation. 

Knowing (7.6) we get from (7.4) and (7.5) that lim ^A^ = 0, lim = -5p(0)-£/r(0), 

n—foo n— >c» " 

and 

Jim / (log (!-£)+ J) = / (log (1 - £) + £) ,,,(s) (7.8) 

for zG C\Fe. 
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We next consider the measures pn defined by 

dpn{s) = -^dvEr,{s). 

It will be our goal to prove that 

sup / dpn < +00 (7.9) 

neN J 

and that the sequence {pn) is tight, that is, 

V£>0:3i2>0:VnGN: / dpn{s)<e, (7.10) 

J\s\>R 

see [6, Chapter 6]. We obtain (7.9) and (7.10) from lower bounds on the absolute values of 

the zeros Zj^n of En- 

Let Nn{r) be the number of non-zero roots of En with absolute value < r. From the 
theory of entire functions, we have that 

Nn{r) < log rnax |/„(z)| - log |/„(0)|, (7.11) 

\z\=er 

see [25, Section 2.5], where 

fn{z) = Z-^''-''En{z). (7.12) 

The integral representations (1.4)-(1.6) and (1.7) imply that 

^ , , (-l)"+in! 1 / e^''^'"dw 

(3n)" 2m Jcj^ [w{w^ - 1)J"+^ 

where Cr is a circle centered at the origin with radius -R > 1. Prom (7.13) it immediately 
follows that 

and so 

_ (-l)"+^n!(3n)^"+^ 

(3n + 2)! • ^^-^^^ 

Using £^n(-s) = Pn{z)e~^^^ + Qn(-2) + Rniz)e^^^ and the fact that the polynomials Pn, Qn, 
and of leading coefficients (—1/2)"'+^, 1, and (—1/2)"'+^ respectively, have their zeros in 
a compact set (independent of n), we easily get that |£'„(z)| < e*^"'^' for every n G N and for 
every \z\ > R with R sufficiently large, say R> Ro>l. Then, from (7.11), (7.12) and (7.14) 
we see that there exists a constant C > so that Nn(r) < Cnr ii r > Ro- Hence, for every j 
with \zj^n\ > -Ro) we have j < Nn{\zj^n\) < Cn\zj^n\j so that 

\Zj,n\ > ^ if \Hn\ > Ro- (7-15) 

This is our required estimate on \zj^n\- 
To estimate 



/ dpn = - y2T~ 



j=l :-J,n| 
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we note that there are at most CuRq zeros with absolute value < Rq. So for j > CuRq, we 
have \zj^n\ > Roi a-nd we can use the estimate (7.15). For j < ChRq, we use the fact that the 
zeros can only accumulate on r^, which implies that there is a 5 > so that \zj^n\ > ^ for 
every j and n. Then for every n G N, 



f l'""^"'^ 1 ~ {Cnf CRo C^n 

J ^P-^ n ^ S^^n ^ j2 - <52 + [c-ni^o] ' 

•> .,=1 j=[Cnilo]+l ■' ^ 



and (7.9) follows. 

Next, we estimate for R> Rq, 



L 



mm 



1 ^ 1 1 ^ . / 1 C^n^ 



j = Nn(R) + l ' j = NniR) + l ^ 



i V— - V — - 

- n ^ i?2 + n ^ p - R^ [CnR] ' 



and this shows that, given £ > 0, we can choose R> Rq sufficiently large so that J|3|>^ dpn{s) < 
£ for every n. Hence (7.10) holds and the sequence (/>„) is tight. 

Now let p* be the limit of a weakly convergent subsequence of (/9„), where Pn ^ P weakly 
means that J fdpn J fdp for every bounded continuous function /. Because of (7.9), we 
have that / dp* < 00. Let du*{s) = \s\'^dp*{s). Then (7.8) implies 

J (log (1 - f ) + du*{s) = J (log (1 - + f ) dpEis) (7.16) 

for every z G C \ F^;. Let R > and let and he,r be the restrictions to {\z\ < R} of v* 
and He, respectively. Then we get by taking the real part and rewriting (7.16) that 

J log \z - s|dz/^(s) = J log \z - s\dpE,R{s) + hR{z), 

where is harmonic in {|2:| < i?}. Since F^; has two-dimensional Lebesgue measure 0, the 
unicity theorem [30, Theorem II. 2.1] applies, showing that and pE,R are equal. As R is 
arbitrary, we have that v* = he and so dp* is) = rKidpEif')- Since the sequence (p„) is tight, 
and YTidpEis) is the only possible limit of a weakly convergent subsequence, the full sequence 
{pn) converges weakly to T^dpE{s), see [6, Chapter 6]. Then it also follows that the measures 
VEn converge to in the sense indicated in the theorem. 

This completes the proof of the theorem. □ 



7.3 Proof of Theorem 2.9 

Proof. Away from Tq we have that one of the terms in (2.27) dominates the others, and 
then (2.27) reduces to (2.25). So it suffices to prove that (2.27) holds uniformly for 2; in a 
neighborhood A of a point zq G Fq \ {zi, Z2, 2:3, z^}. So let zq G Fq \ {zi,Z2,Z2, Z4}. We take 
the contours Tp++ and F^++ so that we can choose a disk A centered at zq that lies in the 
region bounded by Fp++, F^++, F[/m and F^m. 
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Prom (1.3) we know that 

Qn{z) = En{z) - Pn(^)e-3"^ - Rn{z)e^'''. (7.17) 



We already know asymptotic expressions for the polynomials P„, i?„ and the remainder 
function that are valid in A, namely (2.23), (2.24) and the last equation in (2.26). Plugging 
these three formulas into (7.17), we deduce 



,n[2gp{z)+gii{z)-e-3\ogz-3z] 



^n[gp{z)+2gR{z)-l-Z\ogz+2,z\ 



l + O 



l + O 



+ 1 



Using (4.12) and (4.13) we are led to (2.27), uniformly for z e A. 

The proof of (2.28) is similar. We use (2.23), (2.24) and the first equation in (2.25). □ 



7.4 Proof of Theorem 2.10 

Proof. Away from Tp one of the terms in (2.29) is dominating the other, and then (2.29) leads 
to (2.23). So it suffices to consider zq G Tp \ {zi,Z2}, and prove that (2.29) holds uniformly 
for z in a disk A centered at zq- We may choose this disk so that it is disjoint from Pi and 
P2, and also so that it is contained in the region bounded by Pp- and Pp+, see Pigure 13. 

We consider z G A. Then we have as before Pn{z) = U2i{z)e^^^^^^^P^^^^^\ but the relation 
between T and U is U{z) = T{z)Vp}{z), or U{z) = T {z)Vp+ {z)Vp++ depending on whether 
z lies in the region bounded by Pp and Pp- (the —side of Pp) or in the region bounded by 
Pp and Pp+ (the +side). This gives 

U2iiz) = T2i{z) ± z-ie-^^e2("+i)'^^(^)r22(z), 

for z on the ±side of Pp. 

Prom the fact that T = SN with S = 1 + 0{^), we get 

p^{z) = e("+i)(fl^'(^)-«) (^N2i{z) ± z-ie2("+i)<^^(^)-^^7V22(z)) (^1 + O 

+ (^Nuiz) ± z-^e2("+i)'^^(^)3^Ari2(z)) O 

+ (^Nsi{z) ± z-^e2('*+i)^^-(^)3-Ar3_2(z)) O . (7.18) 

Since all functions Njk do not vanish in A, and since = —2, we obtain 

(_2)"+ip„(z) = e("+i)»^(^) 

Now we use (4.12) and the formulas in (6.20) for A''2i and N22, and (2.29) follows. 

The proof of the asymptotic formula (2.30) for is similar. □ 





N2i{z) (^1 + j^ijj ±z-ie2("+i)'^^(^)3^iV22(^) (^1 + (^jj 
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7.5 Proof of Theorem 2.11 

Proof. If we unravel all the transformations for 2; G Ai, we find that Y{z) is a product of no 

less than 10 to 14 matrices, the exact number depending on the region where z is. Recall 
that z belongs to one of the four regions f^^{I), f^^{II), f~^{III), and f~^{IV), see Figure 
12. For z G f^^(III), we have the following, where s = {n + 1)^/^/1(2;), 

/e(n+iK/3 \ 

Y{z) = e'^'/^V^ e-2("+iK/3 S{z)N{z) 

\ e(«+i)^/V 

^^i/2^3z/2 0\ /I -1 \ /s^/^ 0\ /Ai(s) -Ai(aj|s) 0^ 

z-V2e-3^/2 o h i Ai'(s) -col Ai' (cols) 

1/ VO 2e-^'^/V \ 1/ \ 

g-iTT/e 0\ / 1 0\ /^-i/2e-3^/2e(n+i)^p(^) g 0^ 

e'""/^ ° ^ ° ° ^l/2g3^/2g-(n+l)¥.p(^) 

1/ V 1/ V 1, 

1 o\ /I 0^ 

^-lg-3zg2(n+l)<^p(2) 10 1 

1/ VO 2;e-3^e-2("+i)¥'«(^) 1^ 

,g-(n+lK/3 Q \ /e("+i)9^'(^) \ 

Q g2(n+l)<?/3 I I e-("+^)[s^(^)+»«(^)l 

g-(n+i)£/3 y y e*^""'"-^^^-^*^^) y 

(7.19) 

Recall that Pn{z) is the (2, 1) entry of Y{z). We multiply (7.19) with the first unit vector 
(1 0) , and we obtain after some simple calculations, where * denotes an unspecified 
unimportant entry 



* ^ /^e(n+i)(.p(.H^p(.)) /I 0\ /I -1 0\ / .V^Ai(.) 

* / ^ ^ \o */ Vo */ V 



(7.20) 

Similar calculations for z in the other regions near zi lead to exactly the same expression 
(7.20). Then further calculations show that for z G Ai, 

(_2)"+ip„(z) = ^e("+^)(^^^^)+'^^(^» 
r 3 

E ^2fc(^) (Nkiiz) + iz-^e-^'Nk2{z)) (n + lf'^h{zf'^ Ai((n + l)2/3/i(2)) 
.fc=i 

3 

+ E '^2ik(^) (-iVfei(z) + iz-^e-^'Nk2{z)) {n + l)-^/Vi(^)-'/' Ai'((n + 1)^/3/1 (^)) 
ifc=i 

(7.21) 

From the jump condition (6.1) for TV on Fp it easily follows that for each k, 

{Nki{z) + iz-^e-^'Nk2{z))^ = i {Nki{z)+iz-^e-'''Nk2{z))_ , z^Tp. 
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The fourth root in is defined with a cut along Tp, and on Fp there is a jump 

Thus the products 

arc analytic across Tp. Since the functions Nf^j have a fourth root singularity at zi, and 
has a fourth root zero at zi, these functions are non-zero at zi, and therefore also in 
a neighborhood of zi. By shrinking Ai, if necessary, we may assume that these functions are 
without zeros in Ai. Similarly, we have that 

{-Nki{z) + iz-'e-''Nk2iz)) 

is analytic and without zeros in Ai. 

Now we recall that S{z) = I + O (^). So it follows that in the two sums in (7.21) the 
terms with k = 2 are dominating, and we get 

n^/' {N2i{z) + iz-^e~^'N22{z)) h{zf'^ M{{n + 1)^/^1 (^)) + O 
+^-1/6 (-iV2i(z) + zz-ie-=^^iV22(z)) /i(z)-V4Ai'((n + 1)2/^1 (^)) [l + O . (7.22) 

This proves (2.34). We have also shown above that the functions hi and h2 are analytic and 
without zeros in Ai. 

For the asymptotic formula for Qn near zi, we proceed in a similar way. Actually, it is 
easier to do the calculations for z in the region /f^(J/) or /f ^(I). The result is that for 
z e Ai, 

'\ 0\ /I -1 0\ / si/4w|Ai((j2s) 
z-^e-^^ 0|(z % 1 I s-V4a;3Ai'(a;2s) I , (7.23) 




,0 */ \0 */ \ 



where again s = (n + 1)^/^/1(2:). We make use of (4.12) to replace the exponential factor on 
the right-hand side of (7.23) by e("+i)(fF(2)+</'p(2)-32-0. Then, we obtain (2.35) from (7.23) 
in the same way as we obtained (2.34) from (7.20). 

For the asymptotic formula for £"„ near z\ , we may again proceed in a similar way, doing 
the calculations for z in the region f^^{III) or f^^{IV). □ 

7.6 Proof of Corollary 2.12 

Proof. The behavior of the extreme zeros of Pn near zi follows from the asymptotic formula 
(2.34). Indeed, we consider the function 

(_2)n+lg-(n+l)(gpiz)+<pp{z)) 

Fn{t) = ^ ^ w,, , , Pniz), where z = zi + tn-^'^. 
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Then F„ has zeros t^^n = (^iTn — ^i)'^^^^) i/ = 1,2, . . ., and these zeros are ordered by increasing 
absolute value. Because of (2.34) we have that 

F„(t) = Ai ((n + if'^Mz, + (n-^/")) (l + O (i)) 

Since /i is an analytic function with a simple zero at zi and f[{zi) = ci (see (2.32), (2.33)), and 
since hi and ^2 are nonzero analytic functions near zi, we get by expanding these functions 
near zi that 

Fn{t) = Ai (tci + O {n-^"')) + ""'^'^^ Ai' (tei + O (n-'/'^)) + O (n-^'^) . (7.25) 

Now, expanding the Airy function Ai near tci, and observing that the second term of the 
sum in the right-hand side of (7.25) is of order ra~^/^, we get simply 

Fnit) = Ai(tei) + O . (7.26) 

The O-tcrm holds uniformly on compact subsets of the complex f-plane. Prom Hurwitz' 
theorem it follows that for every fixed G N, we have 

lim U n = - — ■ 

n— »oo ci 

Using the fact that — is a simple zero of the Airy function, we obtain from (7.26) that 

This proves (2.40), since z^f^^ = zi + t,y^„n^^/'^. 

The formulas (2.41) and (2.42) for the extreme zeros of Qn and En near zi are obtained 
in a similar way from the asymptotics of Qn and En near zi, given by (2.35) and (2.36), 
respectively. □ 

7.7 Proof of Theorem 2.13 

Proof. We consider z in a compact subset K of Dr. By the third equation in (2.25), we have 



lim -log\Qliz)\ = 2RegQiz), 

n— »oo n 



and by (2.23) and (2.24), 



lim - log \Pn{z)Rn{z)\ = Re {gp{z) + gR{z) - 2i). 

n— ►oo n 

In Dr we may use (4.26), so that 

2Re5Q(z) - Re {gp{z) + gR{z) - 21) = Re {gp^z) - gp{z) + 6z). 
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By (4.12) and (4.13), this is 2Re{(pp - (pR) which is positive in Dr. Thus, APnRn/Q^ is 
exponentially small, and there is a choice of the square root so that 

^Ql{z) - APn{z)Rn{z) = Qn{z) (1 + 0{e-'^)) , 

as n -^^ oo, uniformly for z £ K, for some constant c > 0. Then we choose the — sign in 
(2.44) so that 



^ f . _ -Qn{z) - ^QI{Z) - APn{z)lUz) _ Qn{z) 

2i?„(z) - i?„(z)^' + ^^" 

Now, we use (2.24) and the third equation in (2.25) to obtain 

where the last equality uses (4.26). This proves (2.45) for z G Dr. 

The case z € Dp is similar, but now we have to take the + sign in the formula (2.44) for 
X^. This completes the proof of Theorem 2.13. □ 
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